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ABSTRACT 


An  engineering  model  is  developed  of  the  error  contours  encountered  in 
position  fixing  using  synchronous  time  of  arrival  data.  Normal  distribution  is 
emphasized.  Probabilities  of  fixes  occurring  within  circles  and  ellipses  are 
determined  using  the  natural  oblique  coordinates  associated  with  the  measuring 
system.  Comparisons  are  made  of  three -observer  and  four -observer  configurations. 
Examples  are  provided  and  applications  are  discussed. 
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I.  INTRODUCTION 


In  general,  any  event  generating  a  signal  which  is  sharply  defined  in  time,  and 
which  is  propagated  with  known  velocity  to  three  or  more  fixed  observers,  has  a 
spatial  position  which  can  be  calculated.  Conversely,  from  the  reception  at  a  single 
point  (observer)  of  three  or  more  sharply  defined  time  signals  emanating  from  syn¬ 
chronized  fixed  stations,  a  spatial  position  can  be  determined  (i.e. ,  LORAN). 

The  determination  of  the  position  fix  is  accomplished  by  communicating  the 
synchronous  time  of  arrival  (TOA)  data  to  some  common  point  of  intelligence  where 
the  lines  of  position  and  their  intersection  can  be  calculated. 

The  accuracy  of  calculations  is,  in  general,  directly  proportional  to  the  sharpness 
of  the  time  signature  and  the  separation  (in  space)  of  the  observers  or  the  stations. 

There  have  been  numerous  treatments  (references  1,  2,  3,  4,  5,  6,  and  7)  of 
the  physical  systems  which  provide  measurements  for  calculating  a  position  fix.  Most 
of  these  also  discuss  the  mathematics  of  calculation  of  the  fix,  as  well  as  the  various 
error  factors  encountered  and  their  causes.  Analyses  of  the  impact  of  errors  on  the 
spatial  fix  geometry  has,  however,  for  the  most  part  been  treated  as  incidental  to 
these  other  objectives. 

Therefore,  this  effort  was  undertaken  in  response  to  an  apparent  need  for  the 
derivation  of  a. model  of  the  geometric  analysis  of  errors.  The  prime  purpose  was  to 
make  this  model  thorough  enough  to  explain  the  various  factors  and  parameters  and 
provide  derivations  of  all  functions  and  quantities  necessary  for  understanding  and 
using  such  r  model.  In  the  following  description  of  the  method  of  attack,  most  of  the 
analyses  and  results  obtained  did  not  appear  tc  be  available  today,  as  indicated  by 
research  of  the  references  and  many  other  similar  documents. 

Engineering  simplifications  are  used  for  greater  insight,  provided  the  end  results 
are  not  co-  taminated  by  more  than  5%  error.  Following  a  review  of  the  fundamentals, 
a  new  tool  of  error  geometry  resulting  from  a  uniform  or  constant  error  density  is 
furnished  as  a  means  of  making  comparisons  with  other  density  functions.  Calcula¬ 
tions  and  curves  are  provided  for  the  Quasi -Circular  Error  Probability  (QCEP) 
radius  for  all  values  of  probability  in  addition  to  the  usual  probable  error  of  a  -  50%. 
Determination  is  made  of  the  Elliptical  Error  Probability  (EEP)  ellipse  in  terms  of 
the  probability  of  error,  ratio  of  (space  to  time)  gradients  of  the  base  lines,  and  the 
standard  deviation  of  spatial  displacement  per  base  line.  A  comparison  is  then  made 
of  QCEP's  and  EEP’s. 

Since  both  the  QCEP  and  EEP  analysis  is  performed  for  dependent  (three  observers) 
and  independent  (four  observers)  lines  of  position,  a  comparison  shows  that,  contrary 
to  intuition,  the  three  observers  emerge  superior  to  four. 

Derivation  and  application  of  elliptical  transforms  from  a  rectangular  to  an 
oblique  coordinate  system  is  provided.  Probabilistic/geometrical  aspects  of  inter¬ 
preting  single  data  sets  (versus  the  usual  distribution  cf  large  samples)  are  discussed. 
There  is,  in  effect,  a  correction  factor,  which  is  derived,  for  the  true  divergence  of 
the  oblique  coordinate  system  about  the  fix  point.  Examples  are  also  given  of  applying 
this  model  including  a  possible  deep  space  application. 


II.  FUNDAMENTALS 


The  most  fundamental  tool  involved  is  the  locus  of  points  of  a  given  time  difference, 
i.  e. ,  difference  in  distance  from  a  fixed  pair  of  observers.  This  locus,  as  shown  in 
Figure  l,  is  a  hyperbola. 

The  mathematics  involved  in  this  vital  system  component  are  quite  simple.  If 
the  propagation  from  point  P  to  observers  A  and  8  is  in  a  straight  line  at  a  velocity 
V,  then,  for  a  given  hyperbola, 

-  Dg  =  v(t^  -  tg)  =  constant  (1) 

where  t  a  and  1b  are  the  respective  times  of  arrival  (TOA's)  measured  to  the  same 
(synchronous)  reference,  and  the  distances  Da  and  Db  are  also  straight-line  measure¬ 
ments. 

It  is  important  to  note,  for  the  sake  of  generalization,  that  we  are  not  restricted 
to  any  plane  of  action.  Drawing  the  two  hyperbolas  in  the  plane  of  the  paper  wae 
purely  arbitrary.  The  total  loci  of  points  are  the  surfaces  of  revolution  obtained  by 
revolving  the  hyperbolas,  as  shown  in  Figure  2,  around  the  line  A-B  as  an  axis.  If 
we  wish  to  confine  our  interest  to  any  particular  surface,  which  supposedly  contains 
the  position  fix,  we  may  limit  the  loci  to  those  of  the  lines  of  intersection  of  the 
particular  surface  and  of  the  hyperboloid  of  revolution.  The  major  part  erf  this  report 
is  limited  to  a  plane  representation  of  the  earth,  or  space,  as  such  a  particular 
surface. 

If  the  paths  of  propagation  are  not  straight  lines,  but  have  known  radii  of  curva¬ 
ture,  the  loci  surfaces  may  still  be  defined  and  will,  in  general,  be  hyperbolic  in 
nature.  For  example,  in  the  case  of  an  assumed  sphere  such  as  the  earth,  the  radius 
of  curvature  is  the  radius  of  the  earth  and  the  propagation  paths  are  (assumed)  great 
circles.  The  hyperbolic  surfaces  will,  accordingly,  have  to  be  modified.  In  effect, 
they  would  be  generated  from  a  plane  hyperbolic  figure  whose  distances  Da  and  Dg 
represent  great  circle  arc  distances  instead  of  plane  straight  line  distances. 

While  the  development  of  this  concept  might  be  of  considerable  future  benefit, 
it  is  not  an  objective  of  this  report.  Suffice  it  to  say  that  the  change  in  the  difference 
of  path  lengths  is  not  linear  and  is  given  by 

A(Da-  Db)  =  RfGj  -  02)  -  2R(sin  9^2  -  sin  9^2)  (2) 


where 


R  -  radius  of  the  earth 

and  0?  represent  the  great  circle  distances  of  D^  and  Dg 
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or,  perhaps  more  simply,  starting  with  the  original  hyperbolic  lines,  the  constants 
of  difference  would  be  obtained  from 


-  Dg  =  R (f>2  -  0j)  =  constant  (3) 

It  can  be  shown  that  these  modified  hyperbolic  surfaces  will  intersect  the  earth  in 
such  a  manner  as  to  form  closed  curves  which  are  ellipses,  if  four  observers  form 
two  base  lines  which  are  conjugates  of  each  other.  The  extension  of  such  hyperbolas 
into  ellipses  on  a  spherical  surface  has  been  noted  by  Dr.  E.  A.  Lewis  of  AFCRL. 1 

The  various  formulas  which  have  been  developed  for  the  direct  calculation  of 
position  from  two  or  more  pairs  of  time  differences  are  quite  complicated.  None  of 
these  appears  to  offer  the  accuracy  required  for  most  applications.  One  of  the  most 
promising  methods  currently  in  use  to  obtain  accuracy  is:  obtain  an  approximation  of 
position  using  a  relatively  simple  formula,  calculate  or  look  up  the  time  differences 
which  would  occur  if  this  were  the  true  fix,  and  then  compare  these  time  differences 
with  the  data  sets.  The  resulting  so-called  space/time  error  could  be  reduced  to  as 
small  a  value  as  desired  by  moving  the  position  a  small  amount  in  an  indicated  direc¬ 
tion  and  continuing  iterations  indefinitely.  There  is  usually,  however,  a  practical 
limit  to  the  number  of  such  iterations  allowable. 

Unfortunately,  the  error  encountered  in  this  calculation  is  not  the  only  or  worst 
error.  It  is  the  TOA  measurement  itself  which  is  more  fundamental. 

What  are  the  error  distributions  and  how  do  they  affect  system  engineering  and 
total  system  accuracy?  For  a  given  timing  error  or  timing  error  distribution, 
calculate  geometrical  areas  which  contain  the  measured  point  with  a  given  probability. 

Figures  2  and  3  relate  the  basic  elements  of  the  geometry  involved.  For  a  first 
order  of  accuracy  and  simplification  of  the  model,  it  is  assumed  that  we  are  dealing 
only  with  the  plane  surface  representation  of  the  total  hyperbolic  loci  and  that  the 
propagations  are  straight  lines  in  that  plane. 

The  observers  are  A,  B,  C  and  A,  B,  C,  D  for  the  three-  and  four -observer 
configurations,  respectively.  Also  the  lines  AB,  BC  and  AB,  CD,  respectively,  are 
the  base  lines  of  the  three-  and  four-observer  configurations. 

The  dashed  curves  represent  the  hyperbolas  of  a  constant  difference  of  time  (or 
distance)  as  measured  by  these  observers,  and  selected  to  include  the  point  P.  The 


1E.A.  Lewis,  Geometry  and  First-Order  Error  Statistics  for  Three-  and  Four- 
Station  Hyperbolic  Fixes  on  a  Spherical  Earth,  AFCRL-64-461,  June  1964. 
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asjTuptotes  to  these  curves  are  represented  as  the  straight  lines  from  the  center  cf 
each  bf  jrfe  line  through  the  point  P,  thus  introducing  a  small,  negligible  error  in  the 
angles,  0m,  n-  The  exact  construction  of  the  asymptote  is  accomplished  as  follows: 

a.  The  base  intercept  of  the  hyperbola  is  obtained  from 

2a  =  v  (T.  -  T„)  (4) 


b.  A  semicircle  is  drawn  about  each  base  line  as  a  diameter,  i.e. ,  AB,  BC. 

c.  A  vertical  line  is  erected  from  the  base  line  r,t  point  a  to  the  semicircle 
(Figure  2). 

d.  A  line  is  drawn  from  the  base  line  center  to  this  vertical  intercept  of  the 
semicircle,  and  then  extended  out  as  the  asymptote. 

The  concept  of  the  base  line  error,  e,  is  that  the  quantity  tA  -  *3,  will  have 
errors  ±At  associated  with  it,  also  associated  at  the  base  line  is  a  quantity,  c, 
given  by 

=  v  At_ 

~2 

If,  then,  the  asymptotes  rm,  0m  and  rn,  0n  have  associated  with  them  the  errors 
A  0m,  A  0q  which  result  from  errors  em  and  €a,  the  erosshatched  area  would  repre¬ 
sent  the  area  of  uncertainty  for  given  error  limits.  The  equal  division  or  statistical 
centering  of  A  0  about  the  true  asymptotes  is  implied  only  for  the  sake  of  deriving 
a  correction  factor  (see  Appendix  K).  If  A  0m  and  A  0n  represented  error  limits,  the 
erosshatched  area  would  represent  the  area  of  uncertainty  for  those  limits.  In  the 
sections  to  follow,  except  for  one  portion  which  actually  deals  with  such  a  hard 
limited  case,  the  size  of  9,  and  hence  A  0,  is  assumed  to  follow  a  normal  probability 
distribution.  Another  assumption  made  is  that,  for  purposes  of  engineering  accuracy, 
the  erosshatched  area  is  a  parallelogram  bounded  by  four  displaced  hyperbolic  LOP's 
(lines  of  position). 

In  order  to  measure  cr  evaluate  these  LOP  displacements,  the  oblique  coordinate 
system  m,  n  is  established  at  right  angles  to  rm  and  rn,  respectively,  since  it  is 
assumed  the  LOP's  are  parallel  to  the  respective  rm  or  rn. 

Comparing  the  use  of  this  oblique  coordinate  system  to  the  usual  orthogonal, 
rectangular  coordinate  system  is  believed  to  be  unique  in  that  it  appears  to  be  the 
first  time  that  the  complete  analysis  has  been  performed  on  the  natural  coordinates 
of  the  system.  Secondly,  it  offers  a  certain  advantage  in  eliminating  the  necessity 
for  understanding,  defining,  and  computing,  variances,  standard  deviation,  correla¬ 
tion  coefficient,  etc. ,  required  for  the  use  of  rectangular  coordinates.  Tc  be  more 
analytic,  assume  that  we  are  dealing  with  the  rectangular  coordinate  system  x,  y 
such  that  for  a  true  fix  point  and  given  LOP  displacements,  for  example  m  and  n,  we 
would  have, 

y  =  n 

x  =  m  esc  0  +  n  cot  0 
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It  then  becomes  necessary  to  establish  a  correlation  coefficient  between  x  and  y.  First 
establish  the  variance  in  x  and  variance  in  y.  For  the  case  of  in  and  n,  Independent  of 
each  other,  the  situation  is  not  too  bad,  with, 


a 

y 


u 


cscfc  9  o  2  +  cot2  6  a  2 
m  n 


and  correlation  coefficient  px y  given  by 


COV(x,y)  =  E(xy)  -  E(x)  Efr) 

0  0  0  0 
X  y  X  y 


which  turns  out  to  be 


COT  6  a 


rxy 


V 


CSC2  9  o  2  +  COT2  9  O  2 
m  n 
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Having  determined  these  fundamental  rectangular  parameters,  the  next  step  of 
comparison  would  be  to  integrate  the  probability  density  expression 


1 


2ira  o 

x  y 


xy 


l 

2<'V> 


2xyp 


5L 


o  a 
x  y 


over  the  area  of  interest.  Setting  the  exponent  equal  to  a  constant,  such  as  C2, 
results  in  the  area  being  an  ellipse  with  probability 


which  agrees  with  the  results  obtained  using  the  oblique  system. 

At  this  point  perhaps  the  most  significant  difference  is  in  determining  the  size  of 
the  ellipses  for  a  given  probability;  the  most  significant  advantage  of  the  rectangular 
coordinates  is  that  examination  of  the  equation  of  the  ellipse  can  provide  the  orienta¬ 
tion  and  size  of  ellipse.  By  comparison,  using  oblique  coordinates,  a  differential 
max/min  determination  must  be  made  to  obtain  the  true  spatial  description  of  the 
ellipses , 

Referring  to  the  rectangular  parameters,  when  the  quantities  m  and  □  are  not 
independent  (three  observers),  the  determination  of  variance  x  and  the  resulting 
Pxy,  are  considerably  more  complicated,  resulting  in  much  more  complex  expressions 
for  the  ellipses. 

By  comparison,  measurement  and  analysis  in  the  oblique  system  virtually 
eliminate  the  requirement  for  establishing  these  statistical  parameters  in  the  space 
domain.  As  shown  in  this  report,  once  the  correlation  factor  of  errors  in  the  pairea 
time  differences  (Atn  and  Atm)  is  determined,  and  the  density  function,  f(Atn,  Atm), 
is  known  the  conversion  to  the  spatial  representation  follows  without  further  statistical 
analysis. 

Before  comparing  the  difficulties  of  determining  a  circular  area  which  contains 
the  point  with  a  probability  a ,  a  Quasi-Circular  Error  Probability  (QCEP)  is  defined. 

A  circular  area  of  radius,  R,  regardless  of  the  coordinate  system,  oblique  or  other¬ 
wise,  which  contains  the  point  in  question  with  a  probability  a  is  determined.  This 
differs  from  the  classical  definition  of  CEP^  in  two  respects.  First,  values  of  R  for 
all  values  of  a  from  0  to  1  are  obtained,  rather  than  the  special  case  of  a  =  0. 5  only. 
Second,  the  variates,  whether  x,y  coordinates  or  m,n  coordinates,  are  not  equal  nor 
does  R=/ m2+n2  ;  hence,  there  may  or  may  cot  be  a  Rayleigh  function  involved,  so  it 
is  50%  probable  that  R  exceeds  1. 177  of  some  equivalent  0  .  Nonetheless,  a  joint 
probability  function  throughout  a  circle  is  integrated  and  the  results  are  abbreviated 
QCEP. _ 

3 

Reference  Data  for  Radio  Engineers,  International  Telephone  and  Telegraph 
Corporation^  4th  Edition. 
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Contrary  to  calculations  of  the  QCEP  is  the  oblique  system  (Appendices  C  and  D) 
complication  rather  than  simplification  results  from  utilizing  a  rectangular  coordinate 
system. 

In  tae  case  of  four  observers,  the  rectangular  coordinate  system  requires  a 
solution  of 

R 


where 


ii 


fl  I - o 

4°AV1-',xy 


f2  *  Ta  +  « 
*3  -  Ta  -  « 


f4  *  <A-3-)“2 
4A 


20  O  (1-p  ) 

x  yv  xy  ' 


O  O  (1-p  ) 

x  y'  xy 


0=0 

y  n 


O  =  J CSC2  0  O  2  +  COT2  9  Q  2 
x  w  x  m  x  n 


P  - 
Kxy 


COT  9  0 
x  n 


*CSC2  9  a  2  +  COT2  9  a  2 
x  m  x  n 


•  •mar"* 


whereas,  in  the  oblique  system,  the  calculation 


2v/2?a 


r 


where 


n  -R 

sin  9 


(♦[vj  -«{v2])e*P^-HTj 


dn 


V1  = 


V2 


\HZ  a 


-  (n  cot  9  +yR2-n2) 


m 
sin  0 


V?  a 


-  (n  cot  9*  R2-n2) 


m  . 


For  the  case  of  three  observers,  the  ratio  of  complexity  of  rectangular  over 
oblique  systems  is  considerably  greater.  Both  systems  will  add  a  term  to  the  * 
functions,  hut  for  the  rectangular  coordinates  there  is,  in  effect,  a  correlation 
factor  within  a  correlation  factor,  which  modifies  the  parameters  as  follows 


o  =  o 
y  n 


O  =  CSC  .  - 

x  x  m 


“  9  o  2  +  cot2  9  o  2  +  esc  9  cot  9  o  o 


x  n  m 


5  esc  9  o  o  +  cot2  9  o  2 
2  x  m  n  x  n 


xy 


o  V  esc2  9  o  2  +  cot2  9  o  2  +  esc  9  cot  9  o  o 
nf  xm  xn  x  xmn 


whereas,  for  the  oblique  system,  the  functions  are  changed  only  by  the  addition  of  a 
minus  n// 3/2  oQ  term. 

As  shown  in  Appendix  A,  the  angles  which  form  the  asymptotes  to  the  hyperbolas 
are  given  by 

lv^B  '  *(9 


9  =  cos 

n 


-1 


9  =  cos 

m 


-1 


!  2C 
V*A  "  V 


2C 


(5) 


(6) 


For  all  points  beyond  a  distance  of  one  base  line,  these  asymptotes  are  a  very  good 
approximation  to  the  hyperbolas  themselves,  e.g., 


r  and  r  >  2C 
n  m 


(7) 
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If  greater  accuracy  is  required  for  fixes  which  are  closer,  either  more  precise 
formulas  must  be  used  or  solution  by  construction  can  be  performed.  The  latter 
may  be  performed  by  referring  to  a  preconstructed  system  of  hyperbolas,  provided 
the  incremental  steps  in  time  differences  per  hyperbola  are  fine  enough.  Paren  - 
thetically,  and  as  a  philosophic  rationalization,  it  is  stated  that  the  utility  of  this 
model  is  proportional  in  some  manner  to  the  range.  That  is,  it  is  difficult  to 
visualize  applications  where  rn  or  rm  are  less  than  2C.  One  probably  encounters 
such  things  as:  "I  can  see  the  guy  with  my  naked  eye,  "or,  "  ff  the  person  or  the 
event  is  that  close  to  home,  who  needs  a  system  to  measure  its  location?" 

V/hen  #n  and  9m  have  been  determined,  these  and  the  crossing  angle  8*  are 
related  by 


0x  =  iso0  +  en  -  (0  +  em)  (8> 

where  (ft  is  the  angle  between  the  base  lines.  Perhaps  the  simplest  way  to  determine 
the  distances  rn  and  rm  is  by  construction  on  the  plane  representation  of  the  base  line 
system. 

If  direct  calculation  is  preferred,  it  may  be  obtained  as  follows: 

In  Figure  4,  since  the  line  d  joining  the  centers  of  the  fixed  base  lines  is  iixed, 
the  angles  e  and  f  are  also  fixed.  Then  using  the  law  of  sines 


d 

sin  9 
x 


_ m  _ n _ 

sin  (18O-f-0n)  sin  (^-e) 

sin  (l8O-f-0n) 

sin  9 
x 

= 

sin  9  « 

x 


(9) 

(10) 


(11) 


These  relations  are  equally  applicable  to  either  the  case  of  three  observers  or  four 
observers.  Thus  we  have  performed  either  a  calculation  or  graphical  solution,  or 
combination,  for  all  the  geometric  variables  associated  with  a  given  set  of  fix  data 

Regarding  accuracy  in  the  use  of  plane  geometry,  using  the  spherical  surface 
of  the  earth,  either  a  scaled  sphere  or  globe  or  a  plane  map  f  rojection  of  the  earth 
may  be  used  for  the  construction  and  measurement  of  the  geometric  factors.  The 
former  is  apt  to  be  quite  inconvenient.  Concerning  the  latter,  if  the  map  were 
distortionless,  and  straight  lines  drawn  on  the  map  were  the  true  representation  of 
great  circle  segments  of  the  globe,  then  the  base  lines  and  asymptotes  would  be  true, 
and  the  plane  model  would  be  without  error. 
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Figure  4.  Geometry  of  Angles  Invo^°d 

Since  a  Lambert  conformal  conic  projection  contains  a  distance  measurement 
maximum  error  of  about  1%  to  2%  for  the  ranges  of  concern  herein,  and  since  we  are 
for  most  cases  utilizing  the  asymptotes  as  LOP's  (lines  of  position),  the  errors 
should  be  negligible.  This  is  particularly  true  if  the  crossing  angle,  8*,  is  calculated 
directly  in  terms  of  Lat.  -Long,  and  spherical  angles.  An  enabling  factor  is  the  axiom 
that  the  tangent  to  the  hyperbola  at  the  point  P  bisects  the  angle  subtended  by  the  base 
line,  such  as  APB.  Thus,  if  we  are  seeking  a  statistical  analysis  for  some  assumed 
point,  P,  i»  is  relatively  simple  to  establish  the  angle,  and  62  (see  Figure  3)  of 
the  true  hyperbolas  and  then  obtain 

ex  ■  ei  *  e2  +  eBD 
where,  for  three  observers,  8qd  =  0 
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The  details  acd  derivation  of  such  quantities  as  n,  m,  0q,  8m,  R  and  F are 
given  in  Appendix  A.  As  seen  in  Figure  A-l,  determination  of  these  quantities  is 
independent  of  the  number  of  observers  involved  since  we  are  dealing  with  only  one 
given  pair  of  readings  at  a  time.  In  resume,  the  formulas  for  these  quantities 
follow: 


n  = 


6 

\ 

R' 


r  At  :  m  =  r  At 
n  n  mm 


=  cos 


T  -  T 
-1  V(B  C) 

2C 


sin  8 


1  2  2 

= —  (m  +  n  -  2mn  cos  8  ) 

2  a  X 


(12) 

(13) 

(14) 


x 


where,  in  addition  to  the  quantities  previously  defined,  Ata  is  the  net  timing  error 
made  by  A  and  B  after  subtraction,  and  Atm  is  the  error  of  B  and  C,  or  C  and  D. 


1^  is  a  time  to  space  gradient  whose  value  is  given  by 


r 


n 


r  v 
n 

2C  sin  8 


(15) 


n 

8X  is  the  crossing  angle  formed  by  the  intersecting  hyperbolas  or  their  asymptotes. 
R  is  the  distance  in  the  oblique  axis  system  m,  n  to  any  point  P’  from  the  origin  P. 
C  is  the  base  line  radius  or  one-half  of  the  base  line. 

This  concludes  the  essential  geometry  of  the  fix  system. 


15/16 


in.  DISCUSSION  OF  UNIFORM  PROBABILITY  OF  ERROR  GEOMETRY 

The  probability  of  committing  an  error  of  any  magn  itude  has  not  yet  been 
mentioned.  Of  the  many  possible  distributions  of  errors,  this  discussion  shall  be 
limited  to  two: 


1.  Equal  or  uniform  probability, 

2.  Normal  distribution  of  error. 


The  first  case  will  be  dealt  with  rather  briefly.  While  such  a  distribution  may  never 
occur,  its  analysis  has  academic  value  if  for  no  other  reason  than  to  provide  a 
comparison  or  check  against  results  of  the  second  case  or,  for  that  matter,  any  other 
distribution.  The  details  of  this  analysis  are  given  in  Appendix  B.  The  important 
result  to  be  noted  is  that  for  given  limits  of  error  of  e  the  radius  of  a  circle  which 
would  contain  one-half  of  all  the  measured  fixes,  i.  e. ,  the  QCEP  (quasi-circular 
error  probability)  is  closely  approximated  by  the  following  expressions: 


CASE 

2R  <  AS  <  AS 
n  m 


RADIUS  R 


(16) 


AS  <  2R  <"  AS 
n  m 


AS  <  AS  <  2R 
n  m 


R  = 


R  = 


e/r _ n 

^  4C  sin  0  sin  0 
x  n 


e 

2C 


J _ :*i~ 

^  4  sin^  8  sin* 


a 


(17) 


(18) 


The  QCEP,  for  the  cos  exp  a  =  50%,  means  that  there  i?  a  50%  probability  that 
any  one  particular  calculation  using  data  which  contains  errors  will  produce  a  point  P 
which  lies  within  this  same  circle.  It  can  be  seen  from  Figures  2  and  3,  and  from  the 
previous  discussion,  that  there  is  nothing  inherently  natural  aboi>*  a  circular  area  of 
fix  points.  In  fact,  if  a  continuous  smooth  curve  of  equal  error  probability  were  re¬ 
quired  to  approximate  this  area,  it  would  be  inclined  to  be  an  ellipse.  The  inherent 
benef  .t  of  elliptical  contours  will  be  seen  in  the  second  case  of  normal  distribution 
of  error.  Meanwhile,  we  are  faced  with  a  necessity  of  at  least  being  able  to  compute 
QCEP's  in  order  to  be  compatible  with  other  systems  which  cause  the  event  P  to  occur, 
where  the  probability  of  cause  is  inherently  circular. 

Since  the  second  case  of  a  normally  distributed  probability  of  committing  an 
error  e  is  believed  to  be  more  representative  of  this  system,  it  will  be  treated  in 
considerably  more  detail.  One  example  of  a  component  of  timing  which  may  depart 
from  the  normal  distribution  of  variance  is  the  velocity  of  propagation  of  the  event. 
Actually,  throughout  this  study,  the  velocity  is  assumed  to  be  constant,  i.  e. ,  given 
no  distribution  whatsoever.  However,  if  the  magnitude  of  variance,  with  different 
paths,  is  significant,  it  is  believed  that  its  distribution  is  more  apt  to  be  Rayleigh 
than  normal. 
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IV.  DETERMINATION  OF  STANDARD  DEVIATIONS  OF  LINES  OF  POSITION 
DISPLACEMENTS 


This  section  discusses  the  determination  of  Ob  and  <?rn,  the  standard  deviations 
of  the  statistical  displacements  of  P*  as  measured  along  the  axes,  n,m. 


The  first  step  is  to  establish  the  probability  density  distribution  of  the  timing 
errors  for  a  pair  of  observer  measurements.  This  is  done  in  Appendix  E. 

The  significant  features  are: 


If  o  =  standard  deviation  of  timing  errors  encountered  at  each  observer 

station 

T .  =  exact  true  TOA  at  station  A 
A 


Tfi  =  exact  true  TOA  at  station  B 
The  probability  density  function  of  timing  errors  at  A  and  B  are 

i  r  *a- vi 

f(tA>  =  -  3XP - 2 - 

A  am  L  2c2  J 


=  “ZL  exP 

o/W 


<*B- V 
202 


(19) 


(20) 


and  the  density  function  for  the  pair  is 

1 


f(x)  = 


o  m 

X 


exp 


[-5f| 

u  X  J 


(21) 


where 


x  *A  "  *B  or  *B  "  *C 
M  =  Ta-Tb  or  Tb-Tc 

a  =  /E  a 
x 


This  also  defines  the  important  quantity  <TX  . 


For  convenience  and  reduction  of  terms,  the  quantity  M-X  and  the  quantities  Atn 
and  Atm  as  used  in  Appendix  A  and  elsewhere  are  the  same.  Thus,  to  clarify  and 
standardize  nomenclature,  the  density  functions  for  each  pair  are  written 


f*A  -  tB> 


(22) 
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or 


f(tB  ■  tC)  i  _  1 


f(t 


D  ■  4C' 


f 


(7X/W 


-  exp 


at 


2a 


x 


(23) 


Beferring  again  to  Appendix  A,  and  recalling  that 

v  r 


r  =  n 

n  2C  sin  6 


v  r 

r  =  m 

m  2C  sin  6 


therefore 


m 


n  =  T  At 
n  n 


m  =  F  At 
m  m 


(24) 

(25) 

(26) 
(27) 


^  K  is  Hgain  emphasized  that  this  engineering  simplification  degrades  rapidly  in 
’  accuracy  as  the  distances  rn  and  rm  become  less  than  2C  and/or  as  the  angles  0n 
and  9m  approach  0°  or  180°. 


Consider  the  special  value 


1 


giving  the  special  value 


At  =  a 
n  x 


n„  =  r  a 
°  n  x 


(28) 


(29) 


Using  standard  notation  for  the  normal  density  function  for  variations  in  the  value  of 
n  gives 


f(n)  = 


a  \H?7t 
n 


r  2i 

exp  -JS-g 
20  1 

L  n  J 


(30) 


Substitution  gives 


t(n)  = 


a  J21 t 
n 


exp 


F  2a  2 
n  x 

2a  2 
n 


(31) 
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If  %  is  interpreted  as  that  value  of  n  which  makes  f  (n)  =  0. 242 /<7n  ,  i.  e. ,  the  one 
sigma  value,  then 

0  9  9 

rju^  =  a  *  (32) 

n  x  n 

and  we  have  the  important  relationships 


a  =  r  o 
n  n  x 


(33) 


a  =  r  a 

m  m  x 


(34) 


A  formula  which  gives  a  closer  approximation  for  values  of  rn  or  rm  <  2C  is 
obtained  from  Dr.  E.A.  Lewis10  and  B.  W,  Sitterly  * 


a  = 
m 


JabII—  =  a 
B .  -  B_  x 


2  sin 


B 


2  sin 


VBB 


(33a) 


'A 


where 


"TAB  is  the  standard  deviation  of  the  errors  in'  time -difference  measurement 
made  by  stations  A  and  B,  "  and  is  therefore  equivalent  to  ax  .  Ba  and  Bb  are  the 
''""Respective  bearings  of  the  lines  AP  and  BP  measured  with  respect  to  the  true  north 
vector  at  P.  Thus  Ba  -  Bb  is  the  base  line  angle  APB  subtended  at  P. 


In  fact,  if  we  consider  the  earth  spherical  and  can  measure  BA  and 
Bb  more  readily  and/or  accurately  than  r  and  6  ,  this  computation  of  am  and  an 
where 


a  = 
n 


0  V 
X 


2  sin 


B_.  -  B„ 

t  D 


(34a) 


may  be  preferred  to  the  gradient  rn/m  for 


CTn/m  °x  ^n/m 


used  in  this  study. 


loran,  MIT  Radiation  Laboratory  Series,  Vol  4,  McGraw-Hill,  New  York,  1947. 
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A  more  general  proof  of  the  validity  of  this  concept  follows.  Using  the  probability 
theorem  that  if 


we  have 


v  =  f  (u) 

p(v>  -  pm|$| 


n  =  f(Atn) 

n  =  r  At 
n  n 


(35) 


or 


At  =  #l)  =  w- 
n 


(36) 


hence 


but 


%  -  pHU 


dn  _  r 
d  At  n 


u 


(37) 


so 


and,  if 


p(At  )  =  - -  exp 

a  PE 


[-=?]r- 


p<a  » j  - , — - — 
(an/rp 


exp 


T2At  2 
a  n 


2  a 


n  J 


(Vr^r 


exp 


At 


2a  2/r2 . 

n  n  J 


a 

_  .i 


(38) 


(39) 
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this  becomes  the  identity 


By  symmetry 


P<At  )  » - - - 

n  a  /Sr 


a  = 

x 


exp 


a 

m 

m 


At 


2  <7 


;  2 


x  j 


(40) 


(41) 


Another  relationship  which  introduces  the  important  parameter,  0 ,  is 


Further 


a 

n 

r  a 

n  x 

r 

n 

=  0 

a 

m 

*m  a 

X 

F" 

m 

v  r 

n 

V 

C 

i  ? 

r 

n 

2c  sin  8 
n 

r 

n 

sin  6 

m 

TT 

m 

v  r 

m 

r 

sin  8 

2c  sin  8 


m 


t 


If  g  and  h  are  the  vertical  components  of  rn  and  rm  onto  the  ba$6  lines  and  if 
rn  —  rm  then  for  a  quick  approximation 


(42) 


(43) 


m 


n 

r 

m 


(44) 


This  approximation  is  valid  for  the  geometry  of  both  the  three  and  four  observers,  and 
should  complete  the  basic  model  of  the  system  probability/geometry. 
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V.  DETERMINATION  OF  QUASI-CIRCULAR  ERROR  PROBABILITY  FOR  THE  CASE 

OF  FOUR  OBSERVERS 

From  Figures  2  and  3,  calculate  -he  probability  that  a  given  set  of  time  measure¬ 
ments  will  result  in  a  fix  which  is  displaced  from  its  true  position  P  to  P'.  P*  is 
measured  in  the  oblique  coordinate  system  m,  n  as  the  distance  to  the  vertical 
projections  of  P'  on  n  and  m.  This  is  also  representative  of  the  displacements  of 
the  hyperbolic  LOP's  from  the  true  LOP's. 

At  this  point  the  QCEP  analysis  must  separate  into  two  branches;  as  will  be 
shown,  there  is  significant  difference  between  the  probabilities,  and  hence  radii, 
for  the  case  of  fou*’  versus  three  observers. 

The  case  of  four  observers  is  treated  first  because  of  the  simplification  that  the 
paired  differences  of  time  readings  are  independent  of  each  other,  whereas  such  is 
not  the  case  for  three  observations. 

If  p(n)  and  p(m)  represent  the  respective  probabilities  of  displacements  of  mag¬ 
nitude  n  and  m,  and  if  the  environment  and  readings  of  A  and  B  are  entirely  independ¬ 
ent  of  C  and  D,  then 


f(n,m)  =  f(n)  *f(m) 


(45) 


Also,  for  normal  probability  density  distribution,  we  may  utilize  from  standard 
probability  notation,  the  following: 


f(n)  = 


f  (m) 


a  V2i r 

n 


o  '.'Tn 

m 


exp 


exp 


20 


m 


2a 


m 


or, 


f(n,m)  =  - 


a  a  27r 

n  m 


exp 


n  +  m 


12a  2  2a  2 

n  m  /  j 


The  associated  probabilities  for  incremental  regions  An  and  Am  then  become 

1 


p(n)  = 


a  /27T 
n 


exp 


p(m)  = 


a  v  277 
m 


[-*]* 


m 


(46) 


(47) 


(48) 


(49) 


(50) 
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Consider  now  the  region  A  (which  shall  become  a  circle)  and  all  the  various  combina¬ 
tions  c i  unique  points  n,  m  and  their  probabilities  of  occurrence.  Since  these  are 
mutually  exclusivo  events,  the  total  (collectively  exhaustive)  probability  for  the 
region  A  Is  given  by 


a  =  Yj  f#1)  xf(®)  Am  An 


m  n 


or,  by  well-known  use  of  the  integral  calculus  for  the  continuous  case, 

For  the  special  f^CSP)  case  of  the  region  A  bounded  by  a  circle  of  radius  R 

R  -2  W  r  9  T  r  9  1 


dmdn 


where 


n  r  9  T  r  9  n 

1  r  r  n  in 

a  o~  J  J  ex*>  ~  ~  2  ex? - 2 

m  n  D  .  .  .  2a  2a 

-R  f,(n)  u  n  -»  L  m  J 


i.(p)  =  n  cc«  9  -  wr-E2  sin  6 
1  xv  s 


dmdn 


f2(n)  =  n  cos  9%  +^R2-n2  sin  9  (55) 

A  close  approximation  for  small  values  of  R  sin  9%  was  found  to  be 

r  2  *1  t56) 

exp  -  -^-s-  dn 

L  J 

A  more  exact  expression  developed  in  Appendix  C,  which  is  good  for  all  9  up  to  the 
practical  limit  of  90°,  is  given  by  x 

a  =  — 1 -  f  [*  (x)  -$  (x)]  exp  dx  (57) 

2£/2tTK  h  1  *  L  2K0  J 
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where 


* .  (x)  =  9  |  --  (x  cot  9  +  Jl  -  x^)  I  (58) 

1  (  VTk  ) 

*2(x)  =  (*  cot  9  -7l  -x^)j  (59) 

and  the  symmetry  of  results  is  such  that  solutions  can  be  miide  in  terms  of  that 
variable  whose  standard  deviation  is  the  greater.  The  tabulated  results  follow. 

TABLE  I 


NOTATION  AND  VALDES  OF  £>CEP  VARIABLES 


Figure  5  represents  the  results  of  the  use  of  these  equations  to  determine  the 
circle  R  versus  the  required  probability,  a  .  The  line  drawn  through  the  curves  at 
a  =  50%  gives  the  required  values  of  R,  in  terms  of  On  or  Cm,  for  the  specia1  QCEP 
equivalency  ease. 

Note  that  0,  0n,  and  0*  appear  related  in  some  complex  manner.  If  this  is  true, 
one  of  the  parameters  0  or  6X  could  possibly  be  eliminated  thereby  reducing  the 
number  of  curves  necessary  to  describe  the  system  under  these  conditions.  The 
fact,  however,  that  0  and  0X  are  not  entirely  independent  of  each  other  does  not  in¬ 
validate  their  use  as  parameters.  Values  for  these  curves  were  computed  on  the 
CDC  1604B  computer.  Both  formulas  were  used  for  small  values  of  R  sin  6  with 
very  good  result  comparison.  To  obtain  values  of  the  error  functions  for  $  j  and  $2» 
a  program  was  successfully  written  based  on  the  infinite  series  given  in  Jahnke  and 
Emde  which  resulted  in  good  or  better  accuracy  than  the  values  in  their  table. 


Figure  5(a).  QCEP  Radii  Versus  Probability  a  for  the  Case  of  Four  Observers 
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Figure  5(c).  QCEP  Radii  Versus  Probability  a  for  the  Case  of  Four  Observers 


Figure  5(d).  QCEP  Radii  Versus  Probability  a  for  the  Case  of  Four  Observers 
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VI.  DETERMINATION  OF  THE  QCEP  FOR  THE  CASE  OF  THREE  OBSERVERS 


Compare  the  case  of  four  observers  with  that  of  the  three  observers  In  Figure  2 
where,  in  effect,  observers  Bland  D  are  merged  into  one.  The  situation  with  respect 
to  the  independence  of  pairs  of  time  differences  is  now  changed;  they  are  no  longer 
independent.  In  terms  of  a  coefficient  of  correlation,  or  correlation  factor,  consider 
two  theoretical  extremes  to  help  visualize  the  situation.  First  consider  B  to  be 
errorless,  with  all  of  the  Gaussian  distributed  errors  committed  at  A  and  C.  This 
would  again  represent  an  independent  situation  or  one  of  zero  correlation  factor. 

Next  consider  A  and  C  to  be  errorless,  and  assume  all  of  the  errors  to  be  committed 
by  B.  This  would  represent  a  one-to-one  dependency  with  a  correlation  factor  of 
one.  The  fact  that  the  three  observers  create  errors  with  equal  distribution  makes 
the  actual  correlation  coefficient  of  0. 5  seem  intuitively  feasible.  The  proof  of 
p  =  0. 5  Is  given  in  Appendix  J. 

Also  from  Appendix  J  (Equation  J-32),  the  joint  probability  density  function  in 
the  time  domain  is  given  by 


f<ATN’  ATM>  = 


a  2  ir/3 

x 


exp 


-  2/3 


AT„  +  AT 
N 


M 


-  atm  atn 


(60) 


For  comparison  in  the  space  domain  with  the  case  of  four  observers,  make  the 
following  conversion.  If  incremv  Ms  in  the  space  domain  are  related  to  increments 
in  the  time  domain  by  a  constant,  then 

f(n,m)dAs  =  f(AtQ,  At^)  dAfc 


where  dAs  and  dAj  are  the  respective  two-dimensional  increments.  The  constants 
referred  to  are: 


Then 


n 


r 


n 


dAT 


diTM 


dm 


a 

x 


m 


dA  =  dndm 
s 

dA.t  =  dATNdATM 
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and  by  substitution  into  Equation  (60): 


f(n,  in)  dA 


r  r 

n  m 

a  o  n 

n  m 


exp 


” 

_2i 

(  "2 

2 

+  » 

am  \1 

3I 

,  r 2  o  “ 

r  2o2 

r 

r  a  2  J 

\n  x 

m  x 

n 

m  x  /_ 

dm  dn 

r  r 

m  n 


or 


f(n,m)  = 


o  '  ir/S 
n  m 


exp 


m 


m 


nm 

a  o 

n  m 


(61) 


is  the  probability  density  function  in  the  space  domain,  comparable  to  equation  (48) 
for  the  case  of  four  observers.  Further,  in  comparison,  the  results  of  Appendix  D  are 
as  follows: 


a  = 


+1  r  2  i 

-  /  [  4  (x)  -  4>  (x)]exp  — j 

2/3/27TK  1  2  L  2K02J 


dx 


(62) 


However, 


*2<x) 


sin  6 
_ x_ 

SW2K 

sin  6 
_ x 

/372K 


r(x  cot  6 -VTT7F)  — x . 1 

L  2(3/3/ 2  K  J 


(63) 


(64) 


and  the  symmetry  prevails  such  that  the  conditions  and  values  of  Table  I  apply  directly 
to  this  analysis. 

Figure  6  presents  the  results  of  these  equations  to  determine  the  circle  R  versus 
the  required  probability  a  . 


Referring  to  the  curves  of  both  Figures  5  and  6,  the  abscissae  are  K  and  crg 
where  CTa  represents  om  if  On/om  >  1  and  also  represerts  an  if  ^m/^n  >  1*  provided 
either  ratio  is  equal  to  the  value  of  the  (3  for  the  curves. 


Figure  6(a).  QCEF  .ladii  Versus  Probability  a  for  the  Case  of  Th.  ij©  Observers 
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VII.  COMPARISON  OF  RESULTS  OF  FOUR  OBSERVERS  WITH  THREE  OBSERVERS 

A  comparison  vill  now  be  made  of  three  to  four  observers!  If  the  cor¬ 
responding  £  curves  of  probability  for  four  observers  are  placed  on  top  of  the  curves 
for  three  observers,  there  is  a  nearly  constant  shift  of  the  former  to  the  left  or  to 
greater  radii.  Table  II  gives  a  spectrum  sampling  and  comparison  of  values  through¬ 
out  the  systems  of  curves.  A  fairly  constant  ratio  of  R4/R3  is  obtained  for  a  given  0  . 
The  averages  decrease  with  increased  /3  with  some  indication  of  approaching  1. 0 
as  0— •»  00  .  No  effort  was  made  to  deduce  or  prove  this  theory. 

TABLE  II 

COMPARISON  OF  QCEP  RADII  FOR  THREE  AND  FOUR  OBSERVERS 


QCEP  RADIUS 


RATIO 


Four 

Observers 


Three 

Observers 


R4/B3 


1° 

52. 

38.8 

2° 

26.4 

20.0 

4° 

13.4 

9.75 

8° 

6.86 

5.00 

16° 

3.56 

2.65- 

32° 

1.8 

1.58 

1° 

83.6 

67.8 

2° 

41.8 

33.2 

4° 

21.4 

16.7 

8° 

10.8 

8.50 

16° 

5.48 

4.36 

32° 

2.93 

2.49 

'1*<V  -  IfyfcMi- 


VIII.  COMPARISON  OF  QCEP's  TO  THE  UNIFORM  PROBABILITY  OF  ERROR 

Another  Interesting  comparison  results  with  reference  to  Equations  (16),  (17), 
and  (18).  While  these  referred  to  a  somewhat  academic  situation  of  uniform  or  equal 
probability  of  error  magnitude,  the  following  analysis  indicates  that  these  formulas 
may  have  more  value  than  passing  academic  interest. 

Since  in  Equations  (16),  (17),  and  (18)  the  value  of  e  is  indeterminate,  or  at  least 
somewhat  arbitrary,  a  very  feasible  value  to  assign  would  be  that  derived  from  a 
timing  error  equal  to  the  standard  deviation  of  Gaussian  timing  errors.  Thus, 
assuming 


e  =  v  o 
x 

immediately  forms  a  common  bond  between  the  two  systems  of  measure, 
this  translation  in  Equation  (16) 

v  a  I  r  r 

p  _  x  4/ _ n  m 

C  ▼  27T  sin  V  sin  B  sin  r 
n  m  x 


/  9  v  r  v  r 

/  “  y _ m  __  y _ n  *_ 

'17  sin  &  2C  sin  V  x  2C  sin  9  x 
x  n  n 


To  continue 


'  •;  a  T  a  r  a 

it  sin  o  m  x  n  x 
x 


U  sin  B  ffm  CTn 
!  x 


Also,  using 


fl  =  a 


a  =  &  a 
m  n 


V~?-7r  [R  <  a  <  a  1 

nf  17  sin  o  1  n  m 1 

*  x 


Similarly,  translation  of  Equation  (17)  gives 


2  sin 


[  ff  <  R  <  cr  ] 

n  m 


and  for  Equation  (18) 

R  -  Y\l — T —  +4  K<ffm<R'  W 

2  T2  «ln2  6  “  m 

X 

These  formulas  are  also  adjusted  to  the  condition,  a  requirement  for  consistency,  that, 


Doing  this  makes  it  possible  to  look  up  QCEP's  from  Figures  5  and  6  for  the  same 
&x  and  /3  as  used  in  the  above  equations.  This  was  done  for  the  set  of  points  0x.  0 
given  in  Table  in.  Since  reading  the  curves  of  Figure  5(f)  and  Figure  6(f)  above  32° 
results  in  significant  inaccuracies  and  since  formulas  (16),  (17),  and  (18)  incorporate 
approximations,  the  final  tabulated  results  should  be  viewed  with  a  grain  of  salt. 
However,  looking  at  the  ratio  of  QCEP's  of  three  Gaussian  observers  to  the  hard 
limited  equal  probability  error  makers  (R3g/Ru)  there  is  a  strong  indication  that  for 
rapid  calculation  the  use  of 

R3G  =  1,3  Ru 

would  give  figures  within  10%  error.  Ru  would  have  to  be  calculated  in  accordance 
with  Equations  (68),  (69),  and  (70).  Another  relationship  resulting  from  the  above 
substitutions  is 

AS  =  20 
n  n 

AS  =20 
m  m 

The  criteria  for  the  use  of  Equations  (68),  (69),  or  (70)  have  been  added  in  terms  of 
an  and  Om  as  shown. 
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TABLE  IH 


COMPARISON  OF  QCEP  RADU  FOR  GAUSSIAN  AND 
UNIFORM  ERROR  DISTRIBUTION 


QCEP  RADIUS  (0  = 
n 

1  Mi. ) 

RATIO 

Ru 

Gaussian  R„ 

l  G 

e 

X 

0 

Uniform 

Four 

Observer 

Three 

Observer 

Equation 

R3G//RU 

2° 

1 

14.3  mi. 

27.5 

19.5 

70 

1.36 

8° 

1 

3.7  mi. 

6.8 

5.0 

70 

1.35 

40° 

1 

.  99  mi. 

1.6 

1.38 

68 

1.38 

64° 

1 

.84  mi. 

1.15 

1.14 

68 

1.38 

2° 

3 

43.0  mi. 

57.0 

51.0 

70 

1. 19 

8° 

3 

10.7  mi. 

15.4 

13.0 

70 

1.22 

40° 

3 

2.3  mi. 

3.4 

3.0 

69 

1.28 

64° 

3 

1, 67  mi. 

2.5 

2.32 

69 

1.40 

2° 

5 

71.5  mi. 

95.0 

94.5 

70 

1.32 

8° 

5 

18. 1  mi. 

24.5 

22.3 

70 

1.23 

40° 

5 

3.9  mi. 

5.3 

5.0 

69 

1.28 

64° 

5 

2.78  mi. 

4.0 

3.76 

69 

1.35 
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EX.  AN  EXAMPLE  OF  QCEP  RADIUS  DETERMINATION 

The  following  example  is  provided  to  tie  some  of  these  things  together  in  an 
application.  Assume  an  electromagnetic  phenomenon  as  the  system  medium. 
Further  assumptions: 

1.  The  base  lines  =  200  miles  ■  2C 

2.  lB  =  130°  (Fig.  5) 

3.  £_e  =  25°  (Fig.  5) 

4.  ,if  =  25°  (Fig.  5) 

5.  A's  clock  reads  3.29459025  sec. 

6.  B's  clock  reads  3.2968855  sec. 

7.  C's  clock  reads  3.29337755  sec. 

The  first  calculation  is  to  obtain  the  line,  d,  connecting  base  line  centers 

,  n  sin  <f>  sin  50°  101  . 

d  =  C  — : — -  =  100  — ; — ^37.  -  181,1  miles 
sin  e  sin  2u° 

Next  using  Equations  (5)  and  (6),  the  angles  0n  and  9m  are  determined  as  follows: 

t^  -  t  g  =  901.7  usee, 
tg  -  t  c  =311.0  usee. 


6  =  cos 

n 


v  =  velocity  of  light  =-  0.186><10 

‘b> 


-1  v(tA 


2C 


=  cos 


1  0.  186  x  io6  x  902  X  io'6 


200 


0  =  cos-1  0.8386  =  33° 


6  -  cos 

m 


-1  v(tB  " 
2C 


m 


cos  *  0. 2588 


75 


cos 

o 


-1  0.186  x  io6  x  3ii  x  io"6 


200 


Then,  from  Equation  (8) 


0  r  =  180°  +  33°  -  (130°  +  75°) 

0  =  8° 
x 
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-■ r&z*.  * 


Equations  (10)  and  (1 1)  give 


d  sin  (8  -  e) 

m  ’ 

rn  sin 


r  =  181- 1  sl°.5u?-  -  1000  miles 
D  sin  8° 


r  = 

m 


d  sin  {6  +  f) 
sin  B 


m 


181. 1  sin  58 
sin  8° 


-  1107  miles 


From  Equations  (24)  and  (25) 


v  r 


f  _  _ n 

n  2C  sin~F 


r  =  P- 186X10  XI0QQ  *  1. 70S  *106  miles  per  ...cold 
“  200  sin  33° 

v  r 


r  =  ^  X. 1107  =  1.065  >  106  miles  per  second 

m  200  sin 


As  a  final  supposition,  assume  that  statistical  data  available  supports  a  figure  of 
O-0. 7C7  microsecond  per  observer, 

or 


o 

x 

then 

o 

n 

a 


V"2  <7=1  usee,  per  pair 


T  O  =  1.706  miles 
n  x 

r  a  =  1.065  miles 
m  m  x 
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And  the  last  system  parameter  ft  is 


ft 


l.\  _  1.708 

r~  1.065 
m 


1.603 


If  then  we  wish  to  determine  the  QCEP  for  three  observers  for  a  =  0. 5,  reading 
from  Figure  6  (f }  by  interpolating  between  ft  =  i  and  ft  -  2  we  get 


R  =  vHSOa  =  7.07  0  =  7.07  Xj.  065 

e  m 

R  -  7.  53  miles 

as  the  radius  of  the  circle  which  contains  the  true  point  of  emission  with  a  probability 
of  0.5. 

To  summarize,  if  we  are  given  the  fix  position  computed  from  a  data  set  of  rime 
differences  and  can  calculate  the  system  standard  deviations  (aa,  am)  and  crossing 
angle  ^x*  the  formulas  given  heretofore  or  the  set  curves  of  Figures  5  and  6  can  be 
utilized  to  obtain  the  size  of  the  circle  of  error  for  any  given  probability  of  error. 
Conversely,  if  one  is  speculating  about  a  particular  circle,  the  probability  of  a  true 
fix  being  within  that  circle  can  also  be  read  from  these  curves ,  While  sufficient  data 
has  not  yet  beer,  analyzed  to  establish  practical  limits,  i'  is  believed  that  for 

10°  <  6  and  9  <  170° 

n  m 

r  and  r  >  2C 
n  m 

This  process  will  not  result  in  errors  greater  than  5%  in  the  calculation  of  er/  ors. 
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X,  DETERMINATION  OF  THE  ELLIPTICAL  ERROR  PROBABILITY  SURFACE 
FOR  FOUR  OBSERVERS 


Tc  continue  with  the  analysis  of  error  distribution,  we  will  analyze  the  proba¬ 
bilities  associated  with  areas  bounded  by  an  ellipse.  Since,  as  shown  in  Appendix  F, 
an  ellipse  has  been  found  tor  which  all  points  on  the  ellipse  represent  an  equal 
probability  density  of  fix  error,  this  figure  may  be  very  important  to  systems  applica 
tions.  In  Appendix  F  the  complete  development  of  the  time  domain  to  space  domain 
transformation  is  given.  As  in  the  case  of  the  QCEP,  we  will  deal  first  with  a  four- 
observer  system.  The  significant  results  are: 


1.  Given  the  circle  which  considers  the  probability  of  making  timing  errors 
of  magnitude, 


o  9  o 

AT  ~  +  AT  ^ 
m  n 

whijh  results  in  the  spatial  ellipse 


(71) 


where 


(72) 


by 


2.  The  probability  of  a  fix  being  anywhere  on  or  within  this  ellipse  is  given 

2 


a  = 


27 ICC 
n  m 


J  f  fexP 

J  JA 


m 


9  9 

2a  6  2  a  ‘ 

n  m 


dmdn 


a  ~  1-e 


-A 


(73) 


where 


\HT  a 

X 


(74) 
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For  example,  for  the  "natural  ellipse" 

X  =  1  ,  S2  =  2  0  2 
x 


and 

a  =  1-e"1  =  63.22% 

&  ' 

As  a  more  general  example,  to  find  th<  elliptical  constants  for  a  given  probability 
of  error  for  the  following  case 


0  =1  usee 

x 

=  100  miles 

r  =  1000  miles 
n 

r  ~1107  miles 
m 

6  =  33° 

n 

e  =  75° 

m 

a  =0.5 

results  in  the  same  rn,  Tm  as  in  the  previous  example,  i.  e, , 

T  =  1.706  * 106 
n 

r  =  1.  065  X  106 
m 

As  shown  in  Appendix  E,  the  general  elliptical  equation  is  given  by 


n 


m 


\2  o  2  r 2  2  <t  2  r  2 

x  n  x  m 


In  terms  of  the  major,  minor  axis  concept,  if 


=  1 


a2  =  2X2  a  2  T  2 
x  n 

>  2  0^2  „  2  r  2 
b  =  2X  <7  1 

x  m 


|75) 
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“T  f 


1 


giving 


The  ellipse  would  be 


x  =Vl0^Trb  =  °-834 


a  ■-  0. 834  x  i.  414  x  io"6  x  l.  706  x  io6  =  2.  012  miles 
b  =  0. 834  x  l.  414  x  IO-6  x  i.  065  x  106  =  1. 256  miles 


and  the  ellipse  becomes 


2  2 

_ 5 _  +__E_  =  i 

(2.C12)2  (1.25S)2 

These  are  the  parameters  of  a  conformal  ellipse  drawn  to  orthogonal  axes.  To  see 
the  true  spatial  ellipse  however,  either  a  graphic  projection  onto  the  oblique  axes  will 
have  to  be  performed  or  mathematical  analysis  using  a,  b,  and  0X  to  obtain  an  a'  ,  b* 
and  axes  shift  angle  \j/  . 

This  section  is  devoted  to  a  complete  analysis  and  discussion  of  the  true  spatial 
ellipse  as  it  exists  on  the  true  axes,  m  and  n.  The  details  of  this  analysis  are  given 
in  Appendix  H  with  significant  results  factored  out  here. 

Figure  7  is  presented  as  a i  example  of  a  graphical  analysis.  To  emphasize  the 
elliptical  appearance,  it  was  decided  to  assume  a 


rather  than  carry  through  the  previous  numerical  examp! °  The  assumed  ellipse 
plotted  on  the  orthogonal  set  n,  z  is 


Graphic  transformation  to  the  oblique  axis  set  n,  m  is  then  performed  as  shown  for 
the  sample  point,  o.  The  values  of  z  are  both  used  to  show  the  two  resulting 
transformed  points,  o'.  The  transformation  is  accomplished  by  projecting  o  onto  z, 
then  equilaterally  from  z  onto  m,  thus  making  it  equal  to  the  z  value,  then  perpendicu¬ 
larly  down  from  m  until  it  intersects  the  n  value  of  o  which  is  also  a  perpendicular 
drop  from  the  n  axis.  While  this  solution  was  no:  performed  with  any  great  accuracy, 
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it  is  intended  to  display  the  significance  and  concept  of  the  true  ellipse  on  the  oblique 
axes,  and  to  provide  approximate  confirmation  to  the  analytic  equations.  Measure¬ 
ments  show  the  transformed  major  and  minor  axes  to  be  approximately 


a'  =  7.3 

b'  =  1.  r> 


with  an  axis  shift  of  approximately 

^  =  60° 


If  we  then  let  x  and  y  represent  the  axes  which  contain  a*  and  b*  ,  we  have  the 
quantitative  relationship 


r> 

4 


(1.  85) 


1 


This  then  represents  a  complete  description  of  a  particular  elliptical  contour  of 
constant  probability.  It  can  be  made  as  accurate  as  the  measurements  of  data 
(At,  T)  and  graphical  solution  will  allow. 


The  inconvenience,  if  not  inaccuracy,  of  this  method  strongly  indicated  that  a 
complete  analytical  solution  would  certainly  be  well  worthwhile.  The  results  of  tne 
analysis  given  in  Appendix  H  are  the  following  set  of  transforms  from  the  orthogonal 
to  the  true  oblique  axes. 


R 


where 


Max/ Min  g^2 


.  ab  a 

±  - - —  COE  V 


/ 


X  = 


Thus  when  a  >  b 


1-bVa 2  _  a/b-b/a 

lib/ a  cos  0  2' cos  & 


a  -  vlx  u  T 

x  n 

b  =  VT>  a  r 
x  m 


(76) 


R 


Max 


dr,  *7=! 


(77) 
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Min 


and  when  b  >  a 


1  \  2  ab  co  8  6  ) 

d'Tzn 


p2 


M“  "  2  Bta2  6 


C08J?  I 
2  ! 


(79) 


R 


Mln  j^?e 


1-1*  1F5)  *i  -iR?)-^-'! 


(80) 


Association  of  the  (fluidities  RMaXf  RMim  a',  and  b'  with  the  axes  x  and  y  is 
somewhat  arbitrary.  In  general  with  m  and  n  oriented  as  shown  in  Figure  7,  there 
will  be  an  axis  of  transformation  which  lies  in  the  quadrant  defined  by  n  and  its 
orthogonal  LOP.  If  we  call  this  the  x  axis,  then  it  will  be  the  major  axis  regardless 
of  whether  a  was  longer  or  shorter  than  b,  and  the  value  of  Rmax  would  be  applied 
to  tills  axis;  this  x  axis  will  always  lie  on  oi  between  the  LOP's.  Perhaps  the  previous 
formula  would  carry  more  general  significance  if  it  were  worded 


2 

-Jhr 


R,.  "  R,,. 

Max  Min 


(81) 


The  angle  of  shift  ip  of  the  pseudomajor  axis  to  the  true  major  axis  is  given  by 

2t 


V  t-~> 


cos  xp  = 


l+x‘  / 


(82) 


_4_ra2/i+-JL_\+b2  +2abcos_ 

sto  6  L  \  KZ)  \  /W7 


which  can  also  be  written  (for  a  >  b) 


sin2  6ll  ++.L  -  1  \ 

2  , 

cos  \p  =  - 

\  '  l+x2  / 

+ 1 A  -  — ~2  +4  f1  -Ji~- -vV--—2^ 

If  l+x  a2  \  If  1+x2  /  f1+x 2 


(83) 
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and  for  o  >  a 


2  f 
cos 


9 

sin  9 


hFS) 


-iF5) 


\  r  ,  2a/b  cos  6 
+  1  +4/1  -  —  «,-  +  - - 


^  A** 


whereas  in  terms  of  the  major  axis  transform 


2  ‘V'T? 

COB2  l//  -  - -  —  S  <85> 

2  (Major  Axis  Xform) 

Figure  8  presents  a  few  of  the  curves  for  these  functions.  The  curves  are  pre¬ 
sented  mostly  to  get  a  feel  for  the  positioning  of  the  axes  of  the  true  ellipse. 
The  number  of  curves  calculated  and  number  of  points  per  curve  are  not  consid¬ 
ered  adequate  in  this  report  for  good  engineering,  accuracy.  The  equations  are 
accurate,  however. 

As  with  most  all  special  cases  where  certain  results  are  invitingly  ob¬ 
vious,  it  is  interesting  to  note,  as  proven  in  Appendix  H,  that  when  0^  «  90°, 
i.e.,  the  true  ellipse  and  the  original  ellipse  are  one  and  the  same,  solu¬ 
tion  of  the  above  equation  gives  $  »  0°.  This  result  corroborates  that  no 
shift  is  necessary. 

Also  derived  in  Appendix  H  are  transformation  formulas  for  the  major  and  minor 
axes.  If  R'Max  is  considered  the  transformed,  or  major  axis  value,  then  R*Max  can 
be  determined  readily  from 


RMax 


sin  9 


l1+Vr^*?(1TFi?) 


2  COB  9 


where  a  >  b.  This  is  a  major  axis  transformation  formula.  For  b  >  a  we  consider 
b  and  b'  -he  major  axes,  and 


— l—  j  4  f1  - 

2  sin2  9  {  b  \ 


1  A  f  2a/b  cos  9 

1+x '  7i+? 


8  *  90* 


a/b  or  b/a>i 


Figure  8.  Major  Axis  Shift  (Four  Observers) 
Similarly,  the  minor  axis  transformation  is  given  by  a  >  b 


2  sin 


71 1  i ?(1‘VrT5)+(1+V1T«2')' 


or,  for  b  >  a 


^Min* 


2  sin 


l2» 
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Figures  8,  10  present  some  of  these  curves  of  major  and  minor  axes  transformation. 
The  transformed  minor  axis  values  approach  1  as  a  limit,  1-e. ,  the  ellipse  neve 
gets  ary  wider  than  its  original  value.  The  transformed  major  axis  values,  however, 
approach(a  or  b^sin  6X  as  their  limit,  meaning  that  as  the  crossing  angle  vanishes, 
the  true  spatial  error  ellipse  becomes  infinitely  elongated.  The  choice  of  abscissa 
value  wen  purely  arbitrary,  and  whichever  ratio  was  greater  than  one  was  selected. 

It  is  noted  that  for  a/b,  or  $,  greater  than  5,  the  transforms  are  virtually  at  their 
limits.  In  fact,  (a  or  b^sin  6X  is  a  good  approximation  for  £  >  2  . 


2  3  4 

a/b  or  b/o  >  I 


Figure  9.  Minor  Axis  Transforms  (Four  Observers) 
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XI.  DETERMINATION  OF  THE  EEP  SURFACE  FOR  THREE  OBSERVERS 


Derivation  of  the  true  ellipse  transforms  for  three  observers  appears  in 
Appendix  I,  with  the  following  results.  The  general  formula  for  the  major  and 
minor  axis  is 


and  the  transforms  are  obtained  from 


depending  on  the  relative  sizes  of  a  and  b  in  accordance  with  the  following  table: 


FORMULA  for 


If 

XFMax  is 

XFMin  is 

a>  b 

R/a  (87) 

R/b  (88) 

b  >a 

R/b  (88) 

R/a  (37) 

|  and  b  are  given  by 

;2  (2  -  2| 008  6X  ~ 5?) 

3^a2/b2  -  2|  cos  9^ 

2  (2  ~  zb</a  cos  ~  b2/a*) 
3((j2/a2  -  2b/a  cos  0,\2 


The  rotation  of  the  major  axis,  from  Ita  pseudo  to  true  position  is  given  by 


where  in  Equation  (89)  \p  la  measured  from  the  m  axis  and  in  Equation  (90)  \ p  ia 
measured  from  the  n  axle. 

Contrary  tc  the  case  of  four  observers  where  the  resolved  major  axis  always  lies 
between  the  LOP's,  the  three  observers  case  results  in  a  resolved  major  axis  which 
will,  for  the  most  part,  lie  outside  the  region  defined  by  the  The  only 

exception  to  this  is  when 


1  <  0  <  2  cos  9* 

and 

9  <  60° 
x 


This  suggests,  depending  on  the  cost  of  a  fifth  observer,  chat  a  live-observer  system 
might  produce  an  intersection  of  ellipses  such  that  the  area  of  uncertainty  is  eignifi- 
easily  reduced. 


Determination  of  the  values  of  a  and  b  requires  a  derivation  similar  to  the  case 
of  four  observers.  In  Appendix  G  the  relationship  between  the  size  of  the  (concentric) 
ellipse  and  the  probability  of  containment  is  given  by 


a 


(91) 


Since  this  formula  and  Equation  (73)  are  numerically  the  same  for  given  values  of  X 
or  S,  Figure  11  is  presented  m  a  quick  reference  of  related  values. 

The  CDC  16Q4B  computer  was  again  utilized  to  obtain  the  sample  curves  of 
Figures  12  and  13  showing  the  transform  values  as  a  function  of  a/b  or  /3  .  Compari¬ 
son  of  these  figures  with  Figure  1-2  shows  that  for  9  <  60°  the  maxima  of  the  minor 
axis  transforms  ard  the  minima  of  the  major  axis  transforms  coincide  with  the 
maxima  for  the  curves  of  Figure  1-2.  Further,  the  contours  of  the  minor  axis  trans¬ 
forms  follow  the  contones  of  the  quantity 


1 

I 

| 


whichever  contains  the  quantity  a/b  or  b/a  greater  than  unity. 


! 
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Minor  Axis  Transforms  (T’lree  Observers) 


Figure  13.  Major  Axis  Transforms  (Three  Observers) 

As  in  the  case  of  four  observers,  the  curves  In  Figures  12  and  13  are  not  sufficient 
to  provide  general  engineering  calculations,  although  they  are  reasonably  accurate  for 
the  values  of  0X  chosen. 
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XH.  COMPARISON  OF  EEP  SURFACES  FOR  FOUR  VERSUS  THREE  OBSERVERS 

The  subject  comparison  is  not  quite  as  straightforward  as  that  of  QCEP.  To 
equate  the  iv/o  systems,  the  following  srmrcach  was  vaken.  Given  the  acceptable 
uncertainty  or  probability  a  we  have 


X  .  B  -/to,  £ 


(92) 


This  determines  the  time-  or  space-limited  configuration 


AT  2  +  AT  2 
m  n 
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atm2  +  atk  -  atm  atn 
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4  CBS 


3  OBS 
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2  2 
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2  a  2  2  a  2 
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M2  N2 

2  _  2 
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4  OBS 


3/2  <7 


m 


a  3/2  cr  a 
r  m  n 


^A2 

o  0BS  (93) 


where  upper  ease  (M,  N)  is  ust^  to  distinguish  values  from  m,  n  and  which  in  terms 
of  the  nontransformed  ellipses  is, 
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m  +  _n_ 
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MN 
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where 


a„  =  iPi).  a 


n4 


b4  *  m  Vl 


a.  *  -TV 2  S 
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To  recall  the  geometric  significance  of  this  we  have 


Therefore,  if,  as  in  the  case  of  the  QCEP,  it  is  assumed  that  the  systems  are  such 
that  on  an  average  or  median,  or  any  other  basis 


Then 


a 

a 


4 

3 


JT 

TWi 


-  1,15 


(95) 


This  advantage  for  three  observers  is  not  the  final  conclusion  since  tne  final  trans¬ 
formed  axis  values  r..c  the  conclusive  factors.  To  derive  equal  #'s  the  ratio  desired 
is 


Major  Axis  (3  OBS) 
Major  Axis  (4  OBS) 


j  XF  Max3 

1715  XF  Max. 

4 


Minor  Axis  (3  QBSi 
Minor  Axli”(4  OBS) 


j  XF  Min3 

1715  XF  Min, 
4 


(96) 


(97) 


Figure  14  gives  the  results  of  such  a  ratio  comparison  for  two  values  of  6X  '30°  and 
90°).  It  can  be  shown  that  these  curves  do  a  flip-f/op  at  the  critical  angle  of  0-  =  60°, 

It  is  again  claimed  that  in  systems  applications  the  greatest  concern  is  for  6X  <  60° 
as  covering  the  majority  of  practical  cases.  Further,  it  is  the  elongation  of  the 
major  axis  which  is  the  worst  offender  of,  or  produces  the  most  damaging  effect  in, 
finding  something  or  somebody.  Thus,  speaking  of  curves  A  and  B,  while  the  ellipse 
ior  three  observers  is  a  little  fatter  than  for  four,  the  major  axis  shows  a  considerable 
improvement  fo’-  three  observers. 
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Figure  14.  Comparison  of  Ellipses  for  Three  and  Four  Observers 
Returning  to  the  example  used  for  QCEP, 

O  =  1.706  miles 
n 

0=  1. 065  miles 
in 

£  =  1.603 

6  =  8° 
x 

and  for  threo  observers 

R  =  7.  53  miles 
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We  now  with  to  determine  the  ellipse  (or  a  “0,5  to  order  th»t  direct  comparisons 
can  be  made  with  QCEP.  From  Equation  (92) 


8  *  V*l0*ira  ■  °- 


834 


and 


a  *  WiSO  =  V372  0. 834  x  1. 706  -  1.74  miles 
11 

b  *  /372  SO  «  /372  f. 834  x  1.065  =  1.09  miles 

ID 


Then  for  the  major  axis  of  the  true  ellipse 


and  for  the  minor  axis 


A  *  *  XF_  =  1.74  x?  2 
msx  max 


Rmin  *  bXFmin  =  1‘09Xi*14 


giving 


Pmax  =  i2. 5  miles 
Rmin  "  1*24  miles 


For  which,  as  discussed  below 


1.24  <  7.  53  <  12.5 


or 


R  .  <  R  <  R 
min  c  max 


Thus  we  have  a  complete  engineering  description  of  the  time  spatial  elliptical  cohtour 
of  a  constant  error  probability.  The  position  of  the  point  P  in  question  with  respect 
to  the  system  base  line  gives  rise  to  the  on,  am  and  6X  of  the  system,  and  the  desired 
reliability  of  measurements,  a  ,  gives  rise  to  the  pseudo-elliptical  qualities  a  and  b. 
From  these,  the  true  ellipse  of  the  probability  a  ,  the  position  of  its  major  and  minor 
axes,  the  lvalues  of  these  axes,  and,  if  we  wish,  a  complete  sketch  or  plot  of  any 
desired  ellipse  can  be  described. 

There  is  a  somewhat  crude  yet  interesting  link  between  the  QCEP  and  EEP  which 
may  be  valuable  for  making  system  estimates.  Observing  Figure  15,  it  is  apparent 
that  thera  is  a  circle  of  radius  R  and  an  ellipse  of  major  axis  Rmax  and  minor  axis 
Rmin.  which  by  virtue  of  their  respective  areas  of  consideration  will  produce  the 
same  probability  of  fix. 
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Thus  we  have  the  linkage 


R  ,  <  R  <  R 
min  c  max 


(93) 


Since  it  is  believed  that  the  elliptical  constants  are  more  meaningful  and  more  readily 
calculated  than  the  circle,  If  the  EEP  is  known,  the  circle  R  with  the  same  probability 
can  be  estimated  from  the  above  limits.  Actually,  if  one  wished  to  pursue  this  further, 
a  much  closer  (weighted)  estimate  could  be  developed  such  as 


R  -  kR 
c  max 


where 


k  =  f 


max 


min 


Sitterly,  in  "LORAN"  (see  rei.on  p.21)  has  proposed  for  the  QCEP  case,  when 
converted  to  the  notation  of  this  mode,  that 
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Xffl.  CONSIDERATION  OF  UTILIZING  A  SINGLE  SET  OF  DATA 

As  a  closing  observation  on  this  mcdel,  it  was  felt  that  a  converse  situation 
concerning  the  use  of  a  single  data  set  should  be  discussed  briefly. 


Figure  16,  Effect  of  a  Single  Data  Set  and  Ellipse  Correction 

Consider  Figure  16.  We  have  established  a  model  which  enables  us  to  calculate 
that  a  signal  originating  at  the  Point  P  will  be  located  by  the  system  withm  or  on  the 
a  =  X  ellipse  with  a  nrobabillty  of  X.  Or,  succinctly,  it  tus  signal  {or  experiment) 
were  repeated  from  P,  b&v  a  million  times,  then  500, 000  of  the  measured  points 
would  be  within  the  ellipse  a  =  0. 5, 

Suppose  we  are  faced  with  a  converse  situation.  The  system  receives  a  particular 
signal,  and,  with  errors,  locates  it  at  Pj.  If  it  is  important  to  describe  an  area  which 
will  contain  the  true  point  P  with  a  certain  probability,  a  ,  what  shall  we  dc?  Since 
we  presumably  have  no  other  information  available,  the  first  impulse  is  to  treat  Pi 
as  a  true  point  and  construct,  in  the  manner  of  this  study,  the  ellipse  for  a  =  X  as 
indicated  in  Figure  16,  Another  rationalization  for  this  is  that  (a  priori)  it  is  just 
as  likely  that  we  have  committed  errors  in  one  direction  as  in  another,  or  there 
would  appear  to  be  no  bias  error. 
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Tton  suppose  Motto?  on*  of  these  500, 000  points  appears  to  be  at  Pg.  This  may 
be  eonsktered  as  equally  probable  situation  because  the  probability  density  function  on 
the  periphery  at  the  a  *  X  ellipse  is  a  constant.  Again,  we  wish  to  describe  the  area 
which  will  contain  P  with  a  probability  of  a  *=  X,  and  proceed  as  for  Pi  .  Note  that 
tto  ellipse  about  P2  actually  contains  the  true  point  while  the  ellipse  about  Pi  does  not. 
Is  tills  a  paradox? 

One  possible  answer  would  be  to  search  for  all  possible  true  points  P  whose 
ellipses  would  contain  Pi  as  one  of  its  (assume  500, 000)  measured  points,  A  sample 
scattering  at  such  true  points  or  ellipse  centers  is  shown.  The  outer  limit  cf  such 
a  gathering  would  be  the  egg-shaped  curve  (shown  dotted).  This  curve,  it  is  suggested, 
comes  more  nearly  describing  the  (converse)  areas  based  00  one  sample  measure¬ 
ment  (Pi)  which  would  contain  the  true  point  P  with  probability  a  .  Similarly  the 
normal  or  =  X  ellipse  about  P2  would  be  weighted  by  joining  centers  of  ellipses  which 
certain  Pg  as  a  peripheral  point. 

The  mechanics  for  establishing  this  type  of  an  inverse  function  and  performing 
its  solution  are  beyond  the  scope  ol  this  study.  Suffice  It  to  give  one  word  of  caution 
against  proceeding  along  the  lines  of  a  singular,  discrete,  point  solution.  It  must  be 
remembered  that  the  probability  of  getting  precisely  the  point  Pi,  given  P,  is  an 
infinitesimally  small  number  which  becomes  zero  when  the  point  Pi  consumes  zero 
dimensions.  Thus  a  discussion  or  calculation  of  the  conditional  probability,  p  (P-/P), 
appears  to  be  meaningless.  Likewise,  given  a  measurement  Pi,  calculating  the 
probability  p  (P/Pi)  that  f>  particular  P  is  a  true  point  is  equally  meaningless. 

Hence  we  are  apparently  ixbbed  of  an  opportunity  to  analyze  such  expressions  as, 


P  (Pj/P) 

P  (P/Pj) 


P  (P,  ,  P) 
P(P) 

P  (P,  P,) 


P(P,) 


or  to  apply  Bayes  or  any  other  theorem.  Even  if  we  did  not  have  this  handicap,  there 
is  no  a  posteriori  knowledge  of  the  distribution  of  p(P)  or  p(P<)  unless  some  sort  of 
specific  (say,  target)  information  were  available . 


Fortunately  the  modifying  or  egg -shaping  factor  (if  this  be  the  answer  to  such  a 
requirement)  is  a  relatively  small  percentage  of  the  major  axis  of  the  ellipse. 


XIV.  DIVERGENCE  FACTOR  OF  THE  TRUE  OBLIQUE  COORDINATE  SYSTEM 


Another  and  perhaps  more  plausible  explanation  lies  In  the  basic  premise 
or  In  simplification  of  the  m^del.  As  pointed  out  in  Appendix  A  Equation  (A-12), 


m 


v  r  At 
m  m 

2c  sin  9_ 


The  errors  encountered  in  m,  n  are  n^t  (within  limits  prescribed  for  rm,  rn,  9^ 
the  cause  of  this  so-called  paradox.  Further,  the  elliptical  relationship  (for  three 
observers,  Equation  (G-SO),  Appendix  G) 


2  2 

m  +  n  mn  _  g2 

3/2  a  2  3/2  a  2  3/2  a  a 

m  n  n  m 

using  such  calculations  of  m,  n,  is  perfectly  valid.  However,  this  approximation  is 
apparently  overshadowing  another  assumption  about  which  very  little  has  been  said. 
The  grid  work  for  the  coordinates  of  Appendix  I,  Equation  (I»\),  (1-2),  Figure  1-1, 
etc. ,  assumes  uniform  parallel  lines  which  are  directed  by  the  LOP's  or  direction 
of  the  hyperbolas  at  the  assumed  true  point. 

As  Figure  17  shows,  the  coordinate  grid  systems  derived  from  actual  conditions 
of  Atm  and  Atn  do  not  contain  this  ideal  uniformity.  Introduction  of  a  true  or  non¬ 
constant  grid  about  point  P,  i.  e, ,  one  whicv  follows  the  LOP’s  or  hyperbolas  abcut 
P,  causes  a  variation  of  6X  ,  Looking  at  Figure  17,  the  significant  change  in  the 
mechanics  of  finding  R  is  that  the  true  R  is  PP"  whereas  the  model  In  using  9Xt  m,  and 
n  derives  PP\  In  terms  of  formulas  then,  PP"  is 


R2  =  — \ - [  (m1)2  +  (n,)<i  -  2m!n'  cos  9*  ]  (99) 

sin  0'  x 

x 

Determination  of  the  elliptical  transforms  therefore  requires  finding  the  maximum 
and  minimum  of 

(100) 


dR 

dm* 


=  f(m*.  n»,  9'*) 


where 

^  =  </>( m',  n«) 

As  an  (engineering)  alternative  to  this  complex  problem,  it  is  suggested  that  the 
following  procedure  will  give  a  very  good  approximation  of  the  true  major  axis  values. 
It  can  be  shown  (Appendix  K)  for  a  given  true  point  P  and  displaced  LOP's  that  the 


points  P*,  P",  P,  etc. ,  all  lie  on  a  straight  line.  The  correction  factors  of  the 
original  major  axis  ate  given  by  (see  Figure  K-l), 

e 

c/a  =  / 
x 

6 

d/c  =  -^ 
z 

e 

d/a  =  / 

z 

where  c  Is  the  major  axle  value  obtained  from  application  of  the  transforms,  and  a 
and  d  are  as  shown  below. 


Finally,  as  per  Appendix  K,  this  becomes 


c/a  = 


AO  +  a  e 

9  +  __5_ - SI 

x  2 

d 


9  m 


x  r 


m 


d/c  = 


e 

x  r  r 
n  m 


d/a  = 


6  +J°L 

x  r  r 
_  n _ m 


0  -i-- 

x  rn 


m 


m 


where  n  and  m  ar  e  the  required  coordinate  values  for  Rmax.  The  cha  ige  in  the 
n  inor  axis  value  will  not  be  significant. 

This  modification  or  correction  factor  is,  of  course,  applicable  to  the  basic  model 
irrespective  of  the  so-called  converse  problem.,  i.  e. ,  it  is  applicable  to  the  area  about 
the  true  point  P. 
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XV.  DISCUSSION  OF  SYSTEMS  APPLICATION  WITH  A  DEEP  SPACE  FIXING  EXAMPLE 


While  another  book  could  be  written  on  the  application#  of  this  basic 
model,  comments  here  will  be  limited  to  introductory  general  observations. 

With  regard  to  the  number  of  observers  (stations)  to  be  employed,  the  cost  of 
a  fourth  obsei  ver  must  be  weighed  against  the  'improvement  in  crossing  angle  (9X). 
Included  in  the  cost  factor  is  the  problem  of  communicating  data  from  the  four 
observers  over  greater  distances,  and  xhe  inherent  redaction  in  accuracy  of  tour 
versus  three  observers  for  base  lines  of  equal  or  similar  length.  The  F and  hence  <J 
factors  will  not,  as  a  rule,  demonstrate  any  significant  differences  over  a  given  area 
of  concern. 

The  phenomenology  oi  the  transmission  of  the.  signature  of  the  event,  whether 
longitudinal  (seismic)  or  transverse  (light,  sound,  heat,  radio,  etc.)  requires  only 
that  the  velocity  be  known  ami  that  it  be  either  n  constant  or  vary  with  the  path  in  a 
manner  which  is  known. 

One  example  of  a  sound  or  longitudinal  wave  application  is  the  SOFA?.  system 
which  locates  sharp  disturbances,  such  as  the  boilers  exploring  of  a  ship  sinking  in 
the  sea.  Standard  deviation  values  of  o  were  on  the  order  of  1. 0  to  4, 0  seconds 
depending  on  the  ranges  under  consideration  (up  to  2, 500  miles)  with  accuracy  on 
the  order  of  one  square  mile. 

The  time  resolution  or  sharpness  of  the  signature  of  the  event,  as  stated  before 
has  direct  bearing  on  the  accuracy,  e.g. ,  At,  cr ,  and  ffx,  of  the  system,  hi  a  refined 
model  attention  must  be  paid  to  these  quantities  varying  as  a  function  of  range.  The 
state  of  the  art  in.  generating  and/or  measuring  particular  amplitude/time  or 
frequency/ time  characteristics  has  improved  tremendously  over  t*  ?  last  two  to  three 
decades.  In  the  realm  of  electromagnetic  disturbances,  systems  have  been  developed 
v/hich  obtain  a  o  on  the  order  of  one  microsecond. 

The  signature  must  be  recognized  and  distinguished  from  other  (Intentional  or 
otherwise)  similar  signals.  This  challenge  is  certainly  not  peculiar  to  position¬ 
fixing  systems. 

The  observers  must  use  timing  devices  which  are  precisely  synchronized  with 
each  other  or  which  have  Imown  time  displacements  with  respect  to  each  other. 
Deviations  from  such  conditions  are  charged  directly  to  0^  and  At .  With  the  develop¬ 
ment  of  atomic  clocks  with  drifts  of  1  par*  in  lO1^,  it  is  possible  even  without  syn¬ 
chronizing  communications,  to  stay  within  %  of  1.0  psec. 

Communications  must  be  available  to  transmit  timing  data  to  a  common  point  of 
intelligence.  This  must  be  done  reliably  and  fast  enough  to  keep  up  with  the  average 
rate  of  occurrence  of  events. 

•) 

As  an  ultimate  in  application,  this  model  places  us  in  a  position  to  make  direct 
comparisons  with  position-fixing  systems,  such  as  tracking  space  vehicles  oi 
mlssilec.  These  apparently  utilize  a  single  station  radar  D/F-ing  principle  in 


Mxgtmsdrc  with  luultfv  titlcm  correlation  teonikjisss.  ii  ie  the  belief  of  the  iataor 
t£2t  TOA  systems  will  ultimately  replies  some  B/¥  systems  with  an  order  of  magni- 

improvemeit  «*  accuracy.  Toe  Air  Force  today  has  systems  which  show  TOA 
to  be  definitely  n-9>'  *  ■-  D/F. 

As  an  teounpla,  suppose  we  wished  to  reach  oat  into  deep  space  and  measure 
from  earth  the  position  of  a  space  ship  such  as  Mariner  v.  hare  accuracy  is  extremely 
important.  Suppose  feather  that  we  have  implemented  the  following  synchronous 
system  of  four  satellites  where  1  and  2  are  in  the  equa>;  -rial  plane  and  A  and  B  are  in 
a  polar  plans. 


The  four  coaxial  states  of  existence  which  occur  every  six  hours  are 
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If  each  station  is  acting  as  an  analog  relay  with  transit  time  to  earth  known  to  within 
0. 1  microsecond,  the  four  stations  shown  represent  an  orthogonally  oriented  three- 
dimensional  measuring  system.  Geonv-rlcal  factors  for  one  of  these  states  might 
appear  an. 


Assume,  roughly,  a  mjge  of  1  million  miles  and  a  very  conservative  estimate  of 
10  microseconds  for  ox,  if  we  settle  for  an  uncertainty  factor  of  50%,  we  have  the 
following  parameters 


2C  «  23, 000  miles 
O ,  «  10  x  10"6 

Jfc 

a  =  50% 


Obviously,  and  would  depend  upon  the  relative  orientation  of  the  system  to  the 
ship.  Howeves,  angles  of  600  and  113°  are  fair  representations.  To  estimate  the 
error  magnitude  in  the  plane  of  1-A-shlp 

r  -  r  =  *>• 186  *  ^  »  io'r 

n  m  23  x  io3  x  0.365 


and 


cr  ex  a 


10  x  io"6  x  ic-7 


100  miles 


Since  is  a  critical  pa  aroeter,  a  close  estimate  Is  obtained  by 


e  =  ^lOgoiigooiptoMr  ^1M 

X  io6 


6  =  1°30* 

x 


Locking  up  R  for  the  QCEP  on  Figure  6(f)  gives 

R  =  25. 5  Cf  =  25. 5  x  100  =  2, 550  miles 

•J 

The  elliptical  constants  in  tnis  plane  wex-e  determined  in  conjunction  with  the 
transforms  as 

R  major  =  3, 150  miles 
R  minor  =  100  miles 

To  compare  with  the  following  crude  yet  fundamental  D/F  system  assume  antenna 


rLYING 

OBJECT 


deviations  from  the  true  line  of  propagation  are  of  Gaussian  probability  density  and 
that  a  standard  deviation  of  0. 1°  is  representative,  or 

=  o.i° 

From  this  point  on  we  can  readily  equate  the  two  system,'’  to  the  same  model  by 

a  =  ®  =  0.00175  *  106  =  ;,  75  :  103  miles 

n  r  a 

a  *  a  at  =  0.00175  xio6  =  1.75  x  10°  ml1  as 
m  9  m 


with  a  crossing  angle  of 


n  8, 000  00  -  , 

9  =  ■*—&-  =  28  minutes 

x  106 


34 


Since  the  function  of  antenna  A  should  be  completely  Independent  of  antenna  B,  this 
system  is  similar  to  a  four-observer  TOA  system;  as  a  best  estimate  from  the  curves 
in  Figure  5(f)  for  the  plane  of  A-B-Ship,  the  QCEP  is  given  by 

R  =  300  a  =  1,74  x  io3  x  300  =  5.25  x  io5  mile s 

s 


This  figure  seems  unbelievable  at  first  glance,  like  looking  in  front  of  our  noses  for 
the  traveler.  However,  when  we  examine  the  facts  that  8X  is  less  than  0.5  degree 
and  that  the  beam  uncertainty  is  00  =  0, 1°,  a  total  coll  inear  error  of  1. 0  million 
miles  (2R)  has  some  meaning. 

If  a  correction  cf  the  collinear  error  value  of  R  between  the  earth  and  the  ship 
is  effected  in  accordance  with  determination  of  9y  (see  Appendix  K),we  obtain  a 
value  of 

Ad  A0 

V  •,  +  -TS--r--#-*8+0.l#.  1.5° 

*'  -  2W2  -  OH  *  88.000  mile. 


Finally,  if  we  wish  to  correct  the  major  axle  of  th«  TOA  ellipse,  as  per  Appendix  K 


a/c 


C  + 
x 


"3T0 - 2T~ 

n  .  m 

1  +T~ 


0. 0262 


0.02G2  +2  xiig 


0.995 


A  9  A  9 
p  n  m 

x"  2  “2  _  0.0262  -  0.00012 

F  '  0. u262 


d/c  =  1.005 

Thus  we  have  observed  an  improvement  in  accuracy  of  approximately  iwo  to  three 
orders  of  magnitude.  How  much  better  this  approach  would  be  than  today’s  refined 
space  tracking  D/F  systems  which  may  superimpose  radar  range  or  integrated  CW 
type  of  data,  the  author  has  not  had  sufficient  information  to  determine.  There  is 
also  the  Additional  factor  that  atmospheric  refraction  errors  affect  the  D/F  antennas 
but  do  not  affect  the  TOA  system. 
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APPENDIX  A 


MATHEMATICS  OF  BASIC  GEOMETRIC  ELEMENTS 

This  appendix  is  concerned  with  providing  derivations  of  some  of  the  more  minor 
relationships  used  throughout  the  study. 

1.  Discussion  of 


m 


=  r  At 


m  in 


(A-l) 


In  Figure  A-l,  C  is  the  radius  of  the  circle  and  rm  the  distance  OP.  m  is  the  devia¬ 
tion  of  the  m  axis  or  m  coordinate  of  the  displacement  of  the  true  point  P  due  to  an 
error  A  8  .  a  is  the  base  line  intercept  of  the  hyperbola  through  P.  An  error  in  the 
time  difference  of  observations  at  A  and  B  results  in  the  base  line  error  oi  e  . 

m  is  precisely  given  by 


m  =  r  tan  A  9 
m 

For  A  6  <  1°  we  can  write  with  negligible  error 


tan  A0  = 

C 

Then 

iS 

m  -  r  -pr 
m  C 


(A-2) 


(A -3) 


(A-4) 


In  order  to  get  ef  in  terms  of  e  another  engineering  approximation  is  made,  with  very 
small  error,  that 


Then 


and 


hence 


r  = 


9 

"m 


ef  C 

1  IT 


ef 


m 


r  e 


m 

C  sin  9 

m 


(A-5) 


(A— 6) 


(A-7) 


0-8) 
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Also,  for  perfect  time  readings 

2a  =  vfTA  -  Tb> 

Subtracting  Equation  (A-10)  from  Equation  (A-9)  gives 


Again  substituting 


If  we  choose  to  call 


€  =  i  At 
2  m 


r  v  At 
w  _  m  m 

m  2C  sin  & 

m 


r  v 
m 


r  = 

m  2C  sin  8 


(A-10) 


(A“H) 


(A-12) 


(A-13) 


m 


the  final  result  is 


m  =  r  At 
m  m 


(A-14) 


The  accuracy  of  Equation  (A-5)  and  hence  Equation  (A-12)  dissolves  rapidly  for 

V  a<10°  or  0jr/n  >  170°  . 


2.  To  obtain 


e 


m 


cos 


-1  v»A-tB> 
2C 


(A-15) 


As  discussed  in  the  text  of  tills  study,  it  may  be  necessary  to  determine  and  #n  as 
closely  as  possible  to  either  approximate  a  fix  or  analyze  the  error  characteristics  of 
a  potential  fix  point.  In  the  latter  case,  9m  and  %  would  most  likely  be  geographically 
measured.  In  the  case  of  the  former,  we  would  be  working  with  a  given  set  of  time 
measurements.  In  Figure  A-l 


cos  6  =  a/c 

m 


From  the  hyperbolic  algebra 


thus 


2a  -  v(tA  -  tg) 

-1  v(tA  -  t3> 


(A-16) 


(A-17) 


0  =  cos 

III 


2C 


89 


Since  tA  sad  tg  are  subject  to  the  errors  discussed  elsewhere  the  true  picture  is 
given  by 


cos 


-1 


^A-V 


+  Ad 

m 


(A-.3) 


S.  To  show 


R 


s  in  &  e_ 


(m  +  a  -  2mn  oos 


(A-19) 


This  derivation  does  not  necessarily  apply  to  the  general  problem  of  evaluating  a 
point  P  in  an  oblique  axis  system.  It  may  be  considered  a  special  case  of  an  affine 
transformation  in  which  the  LOP  displacements  from  the  true  LOP's  y  and  z  are 
defined  to  be  parallel  to  y  and  z. 


From  Figure  A-2  the  general  equation  for  the  value  of  m  Is 


m  =  m  +  z  sin  6 
o  x 


(A-20) 


Further  it  can  be  seen  that 


and 


m  =  n  cos  9 
o  x 


m  =  n  cos  +  z  sin  9 ^ 


(A -21) 


Solving  for  z 

m  -  n  oos  9 

z  =  "TO - -  <A~22> 

X 


Also  from  Figure  A-2 

R2  =  n2  +  z2  (A-23) 

but  from  Equation  (A-22) 

z2  s  — -p—  (m2-2mn  cos  9  +  n2  cos2  9  ) 
gin2  9  x  x 


ami,  substituting  In  Equation  (A-23) 

R2  =  — r —  (n*  sin2  9  +  m2  +  n2  cos2  6  -2mn  cos  9  )  (A~24) 

sin2  9  x  n  x 

x 
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APPEND*  1  B 


QCEP  CALCULATIONS  FOR  UNIFORM  ERROR  PROBABILITY  DISTRIBUTION 

Tills  is  an  analysis  of  areas  since  the  probability  of  a  fix  occurring  in  a  g.'ven 
area  is  directly  proportional  to  that  area.  Figures  B-l,  B-2,  and  B-3  illustrate  the 
three  possible  cases.  Quantitatively  stated  these  are: 


CASE  I 

2R  <  AS  <  AS 

n  m 

CASE  II 

AS  <  211  <  AS 

n  m 

C4SE  m 

ASm  <  AS  <  2R 
n  in 

The  difficulty  ia  this  analysis  lies  in  formulating  that  portion  of  the  (QCEP)  circle 
which  lies  within  tbs  error  quadrangle. 

Case  I  is  the  simplest  since  all  of  the  circle  is  contained  within  the  quadrangle. 
Assuming  ASq  <  Asm  for  die  are*  of  the  quadrangle 


Ai  ■ 


AS  AS 
n  m 

bOTS 


(B-l) 


For  the  circle 


Then  for  a  QCEP  of  £0% 


----  vlC 


Ky  r"  l/2  A„ 

to  X 


irR 


,  AS  As 

1  r _ m- 

2  "  sin-# 


(B-2) 

(B-3) 


Solving  for  R 


R  = 


i 


'AS  AS 
n  m 

2tt  sin  6 
x 


(B-4) 


Since 

AS  <  AS 
n  —  m 

assume 

AS  =  y  AS 
m  r  1 

* 

where 

r  =1 

R  = 


V 


Zn  sin  0 


AS 


i  JZHL 

"  v  sin  f 


(B-5) 
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But  i  ince  we  have  the  requirement 


R< 


AS 
_ n 

2 


the  res'll*  is  the  requirement 


or 


V’^x  -  1 

db  ^  1/2 


Tc  is  results  in  the  limits 


y  *  1 
9>  3535' 


and 


r<  1.57 
9  =  90° 


Returning  to  Equation  (B-4)  and  from  Appendix  A 


“i  -  riiet  - 


i 


By  substitution 


R 


,JZi 

y2ir  s 


r  e 

A _  y  Q  y 

sin  9'  0  sin^T 

x  n 


r  c 
m 


m 


R-  eV 


r  r 

n  m _ 

27T  sinT?  Sin'S  sin  W 
x  n  m 


For  Case  n 


AS  <  2R<  AS 
n  m 


A2-2RASn  -  T 
AS  AS 

2RA3  =  --  -7— 7T-- 
n  2  sin  “ 


AS 


R  = 


m 


4  sin  9 


x 


(B-8) 


(B-7) 


(B-8) 


(B-9) 


(E-1G) 

(B-ll) 


(B-12) 
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Finally 


R 


r  e 
n 


4C  gin  9  sin  0 


(B-13) 


cass  m 


AS  <  AS 
n  m 


<  2R 


■I 


This  is  of  interest  because  it  embraces  crossing  angles  of  less  then  39°  which  is 
generally  an  area  of  fix  points  with  range  on  the  order  of  1. 8  bane  lines  and  beyond. 

As  observed  in  Figure  B-3,  a  close  approximation  to  the  enclosed  area  is  given 
by  the  area  contained  in  the  rectangle 


AS  by  W4R2  -  AS  2 
n  »  n 


or 


Aj  =  XiJitf-ZS* 


and  since  Ag  -  — 


- - -  AS  As 

AS  V  4RZ  -  AS  2  = 

n  2  sin  0 


Solving  for  R 


R 


-IS 


f  An  A 


OS 


m  ,  ao  2 


4  sin  0 


a .  +  AS 

2C  n 


By  substitution  via  Equation  (B-8) 


R  = 


w 

2C  f 


m 


4  sin2  9  sin2  6  sin2  6 

x  m  e 


(B-14) 


(B-15) 


(B-16) 


(3-17) 


APPENDIX  C 


QCEP  DERIVATIONS  FOR  A  NORMAL  PROBABILITY  DISTRIBUTION  OF 
ERRORS  FOR  FOUR  OBSERVERS 


A  a  stated  in  the  text,  the  probability  of  a  point  P?  lying  within  a  circle  of 
radius  R  is  given  by 


a  = 


2ir  a  u  J 
n  m  a 


f  f  exp 


• 

r  1 

2 

2  ! 

-n 

exp 

m  1 

2 

2  1 

20 

20 

L  n  J 

L  m  J 

dmdn 


(C-l) 


where  A  is  the  eircie  of  radius  R.  0n  and  0m  are  system  standard  deviations  as 
defined  elsewhere . 

The  immediate  problem  in  performing  the  indicated  double  integration  is  that  we 
are  not  dealing  with  orthogonal  axes  in  which  one  variable  can  be  held  constant  at  some 
known  value  while  summation  is  performed  throughout  the  range  of  the  other  variable. 
In  order  to  utilize  this  accepte  j  method  of  performing  double  integration  one  of  the 
variables  must  be  reduced  to  a  component  which  is  orthogonal  to  the  othor  variable. 

In  Figure  C-l  assume  the  circle  R  provides  tne  desired  probability  a  and  let  z 
be  an  axis  at  right  angles  to  n  .  Any  point  p  on  a  line  parallel  to  z  will  have  a 
constant  value  of  n,  such  ae 


a 


=  ”K 


It  v/as  shown  in  Appendix  A  that  the  relationship  between  the  quantities  n,  m,  ard 
9  is  given  by 

X 


m  =  n  cos  9  +  ^  sin  9 
x  x 

If  this  substitution  is  made  for  m  in  Equation  (A~l)  and  n  is  chosen  as  some 
constant  value,  e.g. ,  nK.  then  variations  in  m  become  a  function  of  variations  in  z 
alone,  which  in  turn  satisfies  the  orthogonal  requirement  for  integration,  i.  e, 

m  =  n„  cos  9  +  z  sin  9  (C-2) 

K  x  x  '  ' 


and 


dm  =  dz  sin  9 


(C-S) 


Substituting  in  Equation  (C-l)  gives 
1 


0!  = 


217  0  a  -  .  , 
n  m  A 


f  f 


-  n 


20' 


n 


exp 


-(zl.  cos  9  +  z  sin  9  ) 

X  X 


21 


20 


m 


dndz  sin  9 
x 

(C-4) 
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Figure  C-l.  Elements  of  QCEP  Geometry 


It  can  be  seen  further  from  Figure  C-l  that  the  limits  of  integration  are 

-R  <  nK  <  +  R 

-  /r5^2  <  z  <  +  V^R2-nK2 

Since  we  are  now  performing  the  integration  first  with  respect  to  z  and  then  witn 
respect  to  n,  the  subscript  K  can  be  dropped. 


For  verification  purposes  it  is  convenient  to  rewrite  the  last  term  os 

_  I  -sin2  9( n  cot  B  +  z)2 
exp  1 - 7  Y- - L- 

2(7 

L  *  J 

anc^  r  +  r  r 

_  sin  9  r  r  -n2  -sin2  9  (a  cot  9  +  z)2 

tt  -  217— ' —  J  — T  “P  - 771 -  1)2  *“  f0  -*) 

m  “  .rJe!-»2  2ffn  J  [  20m 

To  prove  the  correctness  of  this  analysis,  it  is  necessary  that  for  an  infinite  circle 
all  possibilities  are  covered  and 


d2  dn  (C-6) 


a  =  1 


But  Equation  (C-6)  becomes 


sin  9  F  7 

a  =  2ir~a  —  J  J  exP 

n  m  —00  _  00 


M 

2 

2a 
L  n 


|~  -sin2  9( n  cot  9  +  z)2 


i  2a 

L  m 


dz  dn  =  1  (C-7) 


k  «C  U  J 

Since  these  two  integrations  are  to  be  performed  independently,  in  any  order, 
and  since  it  is  easily  shown  that 

I  ex*[i}  *  ■ 1 


it  is  also  necessary  that 


sin  8 


VI  v  0  _« 

in 


/  exp  -ate2  8  ft  cot  g  +  *£  dz  =  i 


(C-9) 


To  check,  let 


when 


and  Equation  (C-9)  becomes 


v2  =  (n  cot  9  +  z)2 
d’  =  dz 


Z  =  °s  ,  V  =  « 


aln  9 
VW  a 


1 exp  [ 


sin2  9  1  ,  sin  9 

2°m  mo 

m 


sin- 8 
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Thus  iVie  necessary  condition  is  met. 


Another  way  at  stating  this  proposition,  which  is  more  general  than  Equation  (C-6), 
Is  to  consider  that  every  N*,  m  (rs  shown  in  Figure  C-l)  traverses  the  values  from 
Mi  to  Mg,  and  to  step  from  the  idea  of  an  orthogonal  component  to  one  of  summing 
all  the  probabilities  associated  with  combinations  of  mi  by  n^  as  mj  goes  from  M2  to 
M2,  i.e. 


where 


,  R  m2  f  2  ‘ 

n1 *»  i  i,  “p  fp 

M1  “  =  n^  coe  ©x  -  z  sin 

Mg  »  f2(n)  =  cob  6x  +  z  sin  ^ 


or  (dropping  the  subscript) 


f ,  (n)  =  n  cos  9  -  sin  6 

1  x  x 

f0(n)  =  n  cos  9  +  J~R^  sin  6 

«  X  X 


I  dndm 


(C-10) 


R  f2(n) 


hr  1  X 

n  m  -R  f 


which  is  the  formula  quoted  in  the  text. 


dmdn  (C-ll) 


Since  this  concept  deals  with  the  well-known  Error  Function  or  Fehlerintegral 
for  which  there  is  no  antiderivative,  and  since  from  a  system  aspect  we  may  be 
dealing  with  small  values  of  9%,  it  is  desirable  from  an  engineering  standpoint  to  try 
to  reduce  this  to  a  single  integral  even  if  only  for  limited  conditions.  One  way  of 
performing  such  an  approximation  is  to  consider  the  integration  of  f(m)  first,  e.g. 


VI?  a 


/  GXP 

fl(u> 


(C-12) 


Referring  to  the  normal  probability  curve  of  Figure  C-2,  the  values  of  fi(n)  and 
f2(n)  are  represented  as  Mi  and  M2.  Since  performing  Equation  (C-12)  is  equivalent 
to  finding  the  prescribed  area  under  this  curve,  |  is  equivalent  to  the  shaded  area 
between  Mj  and  M2.  This  area  can  be, approximated  in  various  ways  such  as  per¬ 
forming  a  series  expansion  about  the  point  {M2  +  Mi)/2.  This,however,  tends  to 
complicate  the  situation.  Therefore,  it  was  decided  that  for  small  values  of  R  sin  9X 


Figure  C-2.  A  Normal  Distribution  Interpretation 


a  reasonable  approximation  is  obtained  by  assuming  the  straight  line  between  p{Mj) 
and  p(M2).  The  area  is  then  the  rectangle 

r  2 
-Mo 

(M2  -  Mp  exp  ! 


20 


1 


m 


plus  the  area  of  the  triangle 

\  (M2  -  Mx)  ^exp 
and  Equation  (C-12)  becomes 


-M. 


20_ 


n 


-  exp 


20 


m  J  i 


\ 

) 


£  - 


27T  O 


m 


(M2-Mi)  exp 


20 


n  J 


2  v-*2 — ^ 


/ 

r .  2i 

r  2” 

\, 

[exp 

-mi 

-exp 

i 

w 

s 

i 

) 

20  2 

2ff  2 

)\ 

\ 

\ 

L  m  J 

L  m  J 

Expanding  and  collecting  terms 

1 


47T  0_ 


m 


2  1 ' 


/ 

r  2  i 

/ 

-Mi 

1  exp 

+  exp 

2  0  2 

V 

\ 

L  m  J 

1- 

-M. 


2-|\ 
7tA) 

&  j  f 
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Substituting  Into  Equation  (C-ll)  and  replacing  M 


a 


2t 

n  m 


R 

-R 


exp 

-fjV) 

+  exp 

-f22(n) 

2 

2  0 

2  0 

L  m  J 

L  m  j  ' 

which  Is  the  formula  quoted  in  the  text. 


da 

(C-13) 


The  limiting  value  of  R  sin  8*  is  somewhat  arbitrary.  It  can  be  seen  from 
Figures  C-l  and  C-2  that  the  worst  error  is  committed  ia  the  region  of  m  =  0  for 

here 


Mg  -  Mj  “  2R  sin  0 


if  it  is  desired  to  hold  this  value  to,  e.g. , 


M2-Mi  »  baB 


then 


bcr 

R  sin  6  < 

X  “*  4 


m 


(C-14) 


While  this  engineering  limit  has  not  been  fully  explored  in  this  study,  it  appears 
from  the  computer  runs  that  good  confirmation  with  the  general  formula  is  obtained 
up  to  a  Am  of  at  least  one  sigma  or  b  =  1. 

Returning  to  Equation  (C-ll)  and  a  more  general  evaluation  of  or  which  is  rigorous 
for  all  8X,  we  proceed  by  using  die  formula  (C-6)  as  follows.  It  is  first  desired  to 
evaluate 

,  2 

-sin  6  (n  cot  6  +  z)  j  ^ 

2 


Let 


X  -  f  exp 


2  0 


m 


u2  =  (n  cot  6  +  z)2 
du  =  dz 


For  limits 


when  z  = 

u  - 

MW 

n  cot  6  ~  J  R2-n2 

MW 

n  cot  8  +  ■/r^-a2 
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and  A  becomes 


n  cot  6  +/  R2-n2 
X  =  f  exp 

n  cot  6  -J R2-n2 


.  2  a  2 
sin  0  u 


2  O’ 


m 


du 


which  can  also  be  written 

n  cot  9 

»-  / 

o 


+  ■/  R2-n2 


exp 


Another  substitution 


n  cot  6  -  vR2-n2 

/  «*P 


2  sin  0  „2 

v  =  - g.  u 


sin2  0  u2 


m 


,2fl  2 

sin  6  u 


m 


du 


du 


20 


m 


and  for  limits 


du  = 


dv 


s 


r^Te 


20 


m 


no 

ro  j.. 


u 

V 

a  cot  6  +  /r2  -n2 

( a  cot  9  +  J R2-n2) 

/I  o 

m 

n  cot  0  -  /r  2  -a2 

(n  cot  0  -  /K?) 

^0 

and  we  again  rewrite  X  as 

no 


x  = 


m 


s  in  & 


V1  V2 
J  exp  [  -v2  ]  dv  -  j  exp  [  -v  ]  riv 


o 


o 


(C-15) 
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fit  the  coefficient  is  rewritten  as 


V7  a  -TB  a 

m  _  2  „  m 

“SETT  =  x  T7ET- 


/Be 


X  “  s- 


» 


3HT^ 


v  v 

|  /  exp  [  ~v2]  dv  - 1  j’exp  [  -v2]  dv 


o  o 

By  definition  the  Error  integral  or  Fehlerintagral  is  given  ty 


(C-1S) 


Therefor® 


*£>  «  —  /  exp[-t2]  dt 

*  o 


V5icr 

x  ‘  TiET> 


Finally,  returning  to  Equation  (C-6) 


a  -  — ~ —  /  [*(^1  -  *{r2)]  exp 

2/2?an  -R 


-n 

2a 


dn 


(C-17) 


(C-18) 


(C-19) 


Since  this  is  not  In  a  form  compatible  with  graphical  or  computer  solution,  the 
following  calculations  are  made. 


Recall  that 


V  =  slnj-  t  t  9  +  y  R2  2  \ 

1  «».V  ' 
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then 


r  =  L 

1  \rYK  \ 


x  cot  9  +  v  1-x 


Similarly 


V  =  L  cot  d  -  /l-x2 

2  /TR  \ 


*(v2)  -  * 


#(V.)  =  *  J  (x  cot  6  +  y^F)  =  *10t) 

— -  fx  cot  6  -  /l  xA  =  *  (x) 

VTK  V  •  2 


Using  £  as  previously  defined,  j3  =  0^/0,^,  we  can  algo  write 


2/32R2K 


=  exp 


and  Equation  (C-19)  becomes 


2/5?  a 


r.jd. 

L  2xr 


/  L*lw  -  *2^ fiXp 


n  -1 


2K  jS  J 


(C»21) 


(C-22) 


(C-2S) 


Substitution  of  0C  =  £  am  ^  CTm/a  =  ^  Into  this  coefficient  completes  the 
transformation  and 


- - ;  f  }*.<x)  ~  4*2fc)]  exp[--^— -1  dx 

nTY  i  1  2  L  2Kp  J 


2/3 /5fK  j 


(C-24) 


This  is  the  formula  quoted  in  the  text. 
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A  word  ok  symmetry,  d,  In  the  original  solution  of  a,  the  initial  integration  had 
bees  performed  for  the  variable  n  Instead  of  m,  ±e  result  would  have  been 


i—  /« P  -- j4-  ♦  k  « W^2) 

VW  _i  2a2/R  |/L'a/R\R  f  11  /J 


(C-25) 


If  the  choice  had  bees  mads  to  let 


x  =  m/B. 

k  =  on2/r.2 

*  ‘  V°n 


the  result  would  have  been 


j— •  /«*[-- £1  *M(,.cot9+/r?)l 

2/3/StTg  L  2P  Kj  L/Hc  \  /-» 


*  I  ita_e  (x  Mt  e 

Lm  \ 


(C-26) 


Comparing  this  with  the  results  of  Equations  (C-21),  (22),  (2U),  and  (24),  it  is 
readily  seen  that  if  the  values  of  K  and  0  in  the  first  solution  are  equal  to  the  parameters 
K  and  £  tu  the  second  solution,  identicri  values  of  a  will  result.  Thus  we  have  a 
choice,  shown  in  the  table  in  the  text,  which  permits  eliminating  the  plotting  of  any 
0  <  1. 

The  abscissae  are  double  scaled  to  give  a  versus  both  K  and  R  where  R  is  in 
terms  of  whichever  o  is  the  lesser,  i,e. 


For  example,  if 


R  =  <7  /4X 

8 


o  >  a  ,  a  a 
n  m  s  m 


It  is  hoped  that  enough  data  is  presented  in  the  curves  of  Figure  5  to  permit  good 
engineering  estimates  for  any  system  requiring  the  QCEP  method  of  procedure. 
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QCEP  DERIVATIONS  FOR  A  NORMAL  PROBABILITY  DISTRIBUTION  OF 
ERROR  FOR  THREE  OBSERVERS 


It  la  shown  in  Appendix  J  that  the  joint  probability  density  function  for  three 
observers  is 


f(n,m)  =  - - exp 

/Staff 

n  m 


/  2  2  ' 

2  /  n  ,  m  mn 

"5  n  72*rr 

\  a  o  n  m 
L  \  n  m  > 


(D-l) 


The  total  summation  of  the  probabilities  associated  with  each  point  m,  n  ia  the  value 
of  a  given  by 

R  r  /  o  o  \  1 


- — - f  J  exp 

^ffanam  -R 


(3/2)o  2  (3/2)o  2  (3/2)0  o 


dmdn 


(L-2) 


wnere 


m  =  f  (n)  =  r.  cos  8  -  y  R  '-n^  sin  8 

*  1  X  3! 

nu  =  (s) r-  n  cos  5^  +  ■»  R  -s  sin  8^ 


a2=-~B- 

(3/2)0 


(D-3) 


(D-4) 


2  m 

v  =  - 5 

(3/2)0  ^ 

m 

Since  we  can  start  the  process  by  holding  n  constant  at  some  value,  e.g. ,  n^ 

R'  v2 

a  =  - — -  r  f  ^  [  -(u4*  +  v?  -  uv)]  x(3/2)o  o  du  dv 

vit  ffo  '  n  m 

n  m  -R'  v. 


(D-5) 


a  *  I  exp  [  -  (U2  +  v2  -uv)  j  dudv 


(D-6) 


/S72  f? 


Therefore,  holding  u  constant 

_  R»  v2  1 

a  “  5?  /  **  /  exp[  -(v2-uv)]  exp  [  -u2]  dv 


(D-7) 


To  evaluate  the  first  integral 


r  z  r  2 

X  =  j  exp  [  -{v  -uv)]  dv 


P-8) 


proceed  by  completing  square 


X  =  J  exp  [ -(v2-uv  +  ~  -  —)]  dv 
V1  4 


P-9) 


=  exp  [  +  -j-  ]  /2exp[-('  -£)]  dv 


y2  =  (v-f ) 


dy  =  dv 

2  y2 

X  =  exp  [  +  ]  J  exp  [  -y2]  dy 

y. 


P- 10) 


Returning  to  a  and  rewriting 


01  TF  f  /  I  exp  t  -y2]  dy  |  exp  [  -  |  U2]  du 

4V7r  -R'  yl  4 


P-11) 
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which  by  definition  of  *  (y)  U 


R1 

a  =  /  I*ty2)  -  *(?!>]  exp[-  Ju2]  du 


p-12) 


Reversing  the  substitution  process 


Rf 

a  =  ?  Vf  /  [*<V  5*  “  «*P  £  -  3/4  u2]  du  P-13) 


^9  “  - ~~ -  m0  =  - — - 0i  cos  0  +  y  R^-n“  sin  6) 

2  V  3/2  a  1  7172  a 


v  s  — 1 — -  0i  co a  0  -  y R2-n2  sin  0 

1  vl72a 

m 

“  ■  ?V?  ft 008  sin  e>  - f  | 

-  *  |  1  * - 0i  cos  0  -  fk^-n1  sin  9)  ~  ~  >  exp[  -0/4)u2]  du 

2)  ( 


And,  a’/ain  for  u, 


ct  ~  t  /  i  *  j" - - - (n  cos  0  +  /r2^2  sin  6} - - — 

-R  |  l/***BV  1  2/B^aJ  ^ 

_# 1 —  (n  cos0  -  Vr2- n2  sin  d) - - - -1 1  exp  j-  I  — - — 1 — d 

L  ^72  am  V  1  27172  an  J  j  [  (3/2)  V  J  ^7: 


7172  a 


AV 

— i—  /  * 

2/2?  a  _r 
n  ^ 


(3/2)a  » 7272  0 

n  J  n 


7172  a  ' 
m 


n  cos  0+7  R  -n  sin  $ 


\  1  2 

"  *  - - - -  (a  C03  0  -  /#? -T?  sin  0) - 5 - exp  f - “r 

7172  a  \  f  27272  a  '  2  0  2 

m  nj  t-  n  J 


(D-16) 


•jSf  jj?  ' 


IWKWH.^XJSV.J.V  -  MV- 


let 


Am  in  the  ears  ct  four  independent  observations,  and  reducing  to  computer  size. 


K  *  J  7B 

m 


x  =  n/R 


0  =  cr  /cr 
n  m 


The  Fehlerlategral  format  becomes 


f  --  R  —  (x  cos  8  ±/l-x2) 
[v  3/2  K  R  '  ' 


sin  0  - 


2/S72"F0 


] 


(D-17) 


1 


2  /IS?  a 


becomes 


exp 


and 


[-*) 


2  0  RvWfC 
becomes  exp 


•  x2  ' 
2K02 . 


dn  =  R  dx 


The  whole  expression  is  then 

+1 
l 

a  ~ 

2/3/5? 


-  « 


SSTT  _j  I  '  /  2  -572K0  J 

I  "" '  1  —  I  x  cos  9  -  ■/  l-x^  sin  0  /  -  ^ — -x  |  exp  j — £—_l 

IV  3/2  K  \  l  2/S72K0J  L  2K02J 


P-18) 

dx 


To  consider  the  symmetry,  had  we  reversed  ice  order  of  solution,  using 

rtl  =  f.  (m)  =  m  cos  6  -  7R2-m2  sin  0 
ll  x  x 

°2  ~  *2  ^  =  m  cos  +  ain  ® 


11C 


WARMS!  S  JBWK* 


and  held  v  constant  while  solving  for  uP  a  direct  reversal  oi' m  for  n  and  c  for  in 
would  have  resulted  as 


K  =  a  2/R 
a 

x  -  m/H 
m  r 

the  final  result  would  be  precisely  the  same  as  JHcpivtiion  (D-18).  Thus  this  solution 
demonstrates  the  same  symmetrical  aspect  as  *  '&  solution  for  four  independent 
observers;  values  of  a  need  only  be  calculated  for  0>  1  . 
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APPENDIX  E 


PROBABILITY  DISTRIBUTION  OF  TIMING  ERRORS  FOR  A  PAIR  OF 

STATION  OBSERVATIONS 

The  problem  is  to  fin.  the  joint  probability  effect  of  a  pair  of  time  oi  servations 
where  each  observation  is  independent  of  the  other  and  each  has  a  normal  distribution 
of  error. 


The  two  stations  Sj  and  S2  are  measuring  the  time  of  arri*  t  the  signature  of 
event,  E.  Assume  the  measurements  are  ti  and  t2  at  Si  and  P  jspectively.  Assume 
the  difference  between  these  two  readings  is 

x  =  t2  -  tj  (E-l) 

The  probability  distributions  can  then  be  expressed  in  terms  of 

p(t2)andp(t2-x) 

We  wish  to  find  the  most  probable  value  of  x  and,  more  generally,  the  probability 
distribution  of  x  with  respect  to  t  .  These  two  events  are  assumed  to  be  entirely  in¬ 
dependent  of  each  other.  Hence  we  are  concerned  with  the  product 

Since  errors  in  measurements  can  be  expressed  as 

tj  “  Tj  +  (E-3) 

'2  =  T2  +  At2  <E-4> 

where  Ti  and  T2  are  the  exact  values  at  Si  and  S2  respectively,  ti  and  (ci-x)  can  take 
on  all  values  within  their  error  spectrum. 

Since  we  are  dealing  with  the  subtraction  of  theae  vime  readings,  what  is  the 
probability  of  subtracting  a  particular  ti  (tu)  from  a  particular  t2  (t2i)  resulting  in 
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The  product  of  the  probabilities  of  occurrence  of  t2i  and  tu  ,  nan  ely, 
p(t2i)  p<til)»  or>  pfct)  *  p(t2i)  *  p(t2i  -  Xj)  Is  the  answer.  But  there  we  many 
possibilities  or  comoinaiioni  or  readings  %2i  and  tu  which  will  give  this  seme  exact 
Xj.  Since  the  occurrence  of  any  one  combination  excludes  the  other b,  the  total  of  all 
possibilities  and  total  probabilities  is  a  summation  process  and  the  total 

P**j)  =  £p  (tg  '  o  ^  -  *j)  (E-5) 


where  i  takes  on  all  possible  values  that  will  result  in  a  particular  time  difference,  Xj. 

Using  a  well theorem  of  probability,  this  csn  be  generalized  for  all  i  by 
the  expression 


PW  -  /  P  (U  F  (t2-x)  dt 


<E-6) 


To  define  the  inr"  ridual  station  probability  densities  let 


P  (to)  =  — - -  exp 

a  -/If 


P  (t. 


r  vvl 

iTJ~\ 

i  r  <i2-x-ti)2i 

-X)  •  — — ;  Kip - j: - 

crV2?  L  2  <r  -J 


(E-7) 

(E~8) 


It  1«  arjumed  there  is  no  reason  why  instrumentation,  human  engineering,  etc. , 
should  result  in  a  differing  from  oue  station  to  the  next. 


Then 


r  1  r  V  T2)2  ’  ^2  "  x  ’  Tl> 

P«  -  /  exp  2  2  1 

-°o  2v  cr  L 


o2 


^2 

(E-9) 


To  evaluate  this  integral  let 


a  =  T„ 


b  =  T  +  x 


Then 


X 

P(X;  =  2"  /  exp 

2jt  <r  -« 


$2 -a)2  -  (tg-b)2 


] 


dt2  (E-10) 
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mtimtcm 


To  simplify  the  exponent. 


2  _2.b2 


/a«b\  ft  + 

It a)  -—3 

“  l/? 


2  a’  ^  f  ft2  +  2ftb  +  b2 


4  0’ 


—~ir  (ft2  +  2ub  +  b  -  t 


7? 


=  -  --V  (a2  -  2ab  +  b2) 
4  c r 


-m2 


and  the  integral  equation  (£-12)  becomes 

o-JW  exp 

then, 


[  /a-b\2l 

[■  \Ta-)  J 


p(x)  =  — x  av? 
2ff  (T 


Replacing  a  and  b  gives, 


PW  = 


2  0-/¥ 


exp 


If  we  now  define  terms  such  that 


I** 


the, refers,  the  most  probable  value  a[  s  ii  x  =  M  fo?  which  p(x)  Is  given  at 

M'T^Tj  (E 

with  a  standard  deviation 

a  =  /Iff 

x 

This  function  is  depicted  In  Figure  E-i  below.  Note  that  while  the  symho.  p&>  is 
used,  the  function  (E-17)  is  a  probability  density  function  for  which  the 
current  trend  Is  to  use  the  notation  f(x). 


Figure  E-l.  Graphical  Results  of  p(x) 
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APPENDIX  F 


<*s?  «■ 


■*$-- 


EEP  (ELLIPTICAL  EBP  OR  PROBABILITY) 

CALCULATIONS  FOR  A  NORMAL  ERROR  PROBABILITY  DISTRIBUTION 

FOR  FOUR  OBSERVERS 

A3  stated  in  the  text,  there  is  nothing  inherently  natural  in  this  type  of  system 
about  a  circular  area  of  probable  occurrence;  and,  as  further  rtated,  If  we  can  break 
free  of  QCEP  compatibility  requirements,  we  can  enjoy  distinct  advantages  in  inherent 
elliptical  contours  of  error. 

To  determine  what  the  inherent  contours  really  are,  we  shall  pursue  the  educated 
guess  that  it  is  possible  to  deter  niDe  the  probability  of  error  directly  from  the  time 
domain. 


As  shown  In  Appendix  E, 

1 


p(x)  = 


a  27r 

X 


exp 


r  (m  -x)2 

[■2a2 

«-  X  -* 


Fcr  the  pair  of  stations  whose  timing  error  would  result  in  deviations  or 
displacements  in  n,  it  is  convenient  to  first  let 


AT  -  M-x 
n 


and 


F^-'J  = 


a  vr27f 

x 


exp 


AT 


2  a 


Likewise,  displacements  in  m  would  result  from 


P(AT  )  = 
r  m 


a  fii 


exp 


AT 


in 


(F-l) 


(F-2) 


(F-3) 


(F~4) 


Recall  from  the  text  that 


a 


p(m)  p(n)  dm  dn 


Since  it  was  also  shown  in  the  text  that 


n  = 

m  ~ 


r  At 
r  At 


(F-5) 


(F— 6) 
F— 7) 
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1 


I 


*****  KXZ&ttCHS 


it  can  be  stated  from  the  principal 


that 


»  <m>  -  T5T~ 

p  (m)  dm  =  p(At  )  dAt  (F-9) 

mm 

p  (n)  dn  =•  p(AtQ)  dAt^  (F-10) 


Hence 


or  =  f  J  p(m)  p{n)  dm  dn  =  j  fT  p(At  )  p{AtQ)  dAtmdAtn 
A 


but  substituting  Equations  (F-3)  and  (F-4) 

r 


a  =  — J—  f  f  exp 

a  2tt  '  t 

x 


<Atn2+%2' 


2  C 


2 


dAt  dAt 
m  n 


(F-ll) 


(F-12) 


Atm  and  Atn  are  independent  of  each  other  and  can  represent  variables  in  the  usual 
orthogonal  relationship  for  purposes  of  geometric  analysis. 

Consider  Figure  F-l  and  the  results  of  letting  the  magnitude  of  timing  errors 
from  0  to  A^  and  0  to  A^  be  represented  by 

At  2  +  At  2  =  S2  (F-13) 

m  n  '  ' 


Figure  F-l.  An  Element  of  EEP  Analysis  in  the  Time  Domain  (Four  Observers) 


\ 
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Thie  vvill  con\ert  Equation  (F - 12)  to 


a 


_i _ 

2V  O  2 
x 


dA 


(F-14) 


where  A  is  bounded  by  the  circle  of  radius  S. 


The  element  of  area  can  be  considered  the  ring  tfor  the  first  Integration) 


whence 


(F-15) 


This  is  the  probability  of  a  fix  lying  within  the  area  resulting  from  the  probability 
of  error  measurements  of  Atj*  and  in  the  time  domain  and  LOP  errors  of  n  and  m 
in  the  space  domain.  To  pursue  this  further  a  is  also  related  to  the  probability  of 
making  LOP  errors  of 


n 


=  r  At 

n  n 


m  - 


r 

m  m 


Recalling  from  Equation  (C-l)  Appendix  C,  *hat 


ot  = 


2r,  a  a 

n  m 


f  f  exp  - 
‘A 


f  2 

a 


2  1 


jl_  +  .. 

2  2  I  ' 

20  “  20  1 

n  m  -J 


drndn 


(F-  16) 
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.  _  - _  .....  - - - - - -  - nwc-  ***!» 


'  {  r*3is*js 


'!?5535^  ■■  ■*-:■■"  **->  ^35.*.  p 


and  that 


o  =  a  T 
n  x  n 

o  =  o  r 


rr< 


x  m 


then 


a  = 


2rr  r  r  a  l 
n  m  x 


I  J  ex P 


dmdn  (F-17) 


since 


2  2 

n  m 

2  r 2  2  r  2 

n  m 


9  9 

At  “  AT  "  c2 
_ m  , _ n_  _  S _ 

2  2  "  2 


(F-18) 


We  have  established  that  taking  the  points  Atn  and  Atm  from  the  circle  S  in  the 
time  domain  results  in  an  elliptical  figure  in  the  space  domain  which  is  characterized 
by  the  quantity 


m 


2  r 


or 


m 


2  r 


m 


2 cr  2  2a  2 
n  m 


a 


It  must  be  remembered  that  a  is  the  probability  of  making  any  combination  of 
timing  errors,  Atn  and  At^,,  up  to  and  including  S.  This  results  in  a  spatial  ellipse 
with  the  same  probability  of  a  fix  being  anywhere  within  or  on  this  ellipse.  The  fact 
that  n  and  m  are  oblique  has  had  nothing  to  do  with  the  analysis  thus  far.  By  letting 


n,.  =  r  s 

Max  n 


m_  _  =  r  s 

Max  m 


(F — 19) 


(F-20) 


the  probability  a  associated  with  any  ellipse 

2 

nMax 


m 


Max 


=  1 


20 


20 


m 


can  be  determined.  However,  in  plotting  this  ellipse,  n  and  m  must  be  the  true  non- 
orthogonal  axes.  Sample  curves  are  presented  in  Figure  F-2. 
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u 


“  >  gk*£k!?l&« «t~-* z%r*m& fir  •  v-f 


Thus  far  we  have  not  placed  any  constraints  or  special  significance  on  the  quantity 
S  itself,  other  than  in  Equation  (F-18).  But  this  represents  a  specific  —  one  might 
say  natural  —  ellipse.  Can  we  make  concentric  expansions  and  contractions  to  this 
ellipse?  To  answer,  it  may  help  to  introduce  the  idea  of  S0  for  this  "natural"  ellipse 
such  that 


S 

o 


■n  a 

X 


(F-?.l) 


representing  that  limit  of  -/At  2  +  At  2  which  is  also  the  limit,  or  boundary,  of  the 
ellipse  m 


m 


20 


m 


1 


To  change  to  a  different  ellipse  of 


n 


m 


r2C„2  *2  20m2 


equivalent  to  the  time  domain  change  of 

S  =  \S 

o 


(F -22) 
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a  -  l-e"’1  =  0. 6322 
resulting  in  the  general  formula 


with  the  spatial  elliptical  relations 


where 


a  =  X  a  =  X  /"S  cr  =  X/2  o  T 
o  n  x  n 


(F  -25) 


(F  -26) 


(F-27) 


b  =  Xb  =  XV2  a  =  X/2  a  T  (F  -28) 

o  m  x  m  ' 

?J b  or  b/a  =  r  JT  or  T  /T  =  j3  [  /3  >  l]  (F-29) 

n  m  m  n  L  J 


Again  note  that  these  are  pseudo-  or  orthogonal -elliptical  quantities  which  must  be 
converted  in  accordance  with  the  discipliner  discussed  in  Appendix  H  to  get  the  true 
spatial  picture. 


124 


-V»  -»*«■  r%  *  •, 


APPENDIX  G 

EEP  CALCULATIONS  FOR  A  NORMAL  ERROR  PROBABILITY  DISTRIBUTION  FOfr 

THREE  OBSERVERS 

To  determine  what  the  inherent  elliptical  contours  of  error  are,  we  shall  pursue 
the  determination  of  error  directly  in  the  time  domain.  For  the  pair  of  stations  whose 
timing  error  would  result  in  deviations  or  displacements  in  n  it  was  shown  in  Appendix 
E  that 


f(x)  = 


a  Jzn 
x 


exp 


(M-x)2 
2  a  2 

L  x  J 


(G— 1) 


Further,  for  this  pair  of  stations  let 

At  =  M-x  =  At„  +  At 


(G-2) 


as  a  simpler  representation  of  the  net  errcr  for  the  time  difference  of  a  particular 
pair  of  stations.  Then 

2 

i  I  01 

f(At  )  =  - —  exp  [  -  — ~  |  (G-3) 


a  V2  TT 
x 


2a 


Likewise,  displacements  in  m  would  result  from  Atm  with  probability  density 


f  (At  )  = 

'  nr 


exp 


At 


m 


2a 


(G-4) 


where 


At 


m 


M-y  =  AtA  +  AtB 


(G-5) 


As  shown  in  Appendix  J,  the  joint  probability  density  function  is  given  by, 

1 


f (At  ,  At  ) 
'  n’  nr 


exp 


VStt  a 


2/3 (At  2  +  At  2  -  At  At 
n  m  n  m 


(G— 6) 


and  the  summation  of  points  within  the  ellipse 


At 


(3/2  XT 


At 


m 


At  At 

_  m  n 

(3/2)0  2  (3/2)0  2 

X  X 


(G-7) 
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gives  a  probability  of  joint  error  for 


°<>tn<A‘N 


in  a  \ 


which  is  given  by 


pfAt^,  At^)  -  a  -  - — -  f  J  exp  - 


r  A  t  2 


•/Jir  a 


t 


At 


m 


At  At 
m  n 


C3/2XT  2  (3/2) a  2  (2/2)0  2 


dAt  dAt 
n  m 


Eo^Uo^G  7hiS  *??/“!  retr  t0  ^  G',1  Whl0h  may  be  ^idered  a  sketch  ot 
Equation  (G-7, .  It  can  be  shown  algebraically  that,  this  is  equivalent  to  a  skewed 

ellipse  onto  axes  Atn  and  Atm'  whose  equation  is  given  by 


(G— 9) 


(At  ')2  (At  ')2 
a  ;  J  1  m  ’ 

r\  '  _  "  » 


a 


The  major  and  minor  axe3  a  and  b  can  be  determined  by 


(G— 10) 


2  2 

J-  =  CQ3  0  sin  ^ 

a2  (3/2)0  2  (3/2 to  9 

x  x 


sin  0  cos  (j> 
(3/2)0  2 


1  __  1  -  sin  cos  0 
a2  (3/2)0  2 


(G— 1 1) 


Imilarly, 


$  is  determined  by 


2 

Si  — 


(3/2)  C 


...  2 


1  -  sin  0  con  $ 

(3/3)0  2 
X 

1  +  sin  $  cos  (J 


$  =  arc  tan 


[3  o  2  3  2\ 

U  x  2  °x  ) 


# 


(G-12) 


O  2  -  |  O  2  |  +  1 
x  2  x  / 


=  45 


arc  tan  1 
o 
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Returning  to  Equation  (G-10),  the  skewed  ellipse  becomes 


(Ay)2 

3a2  c2  '  ' 


(G-13) 


Using  77  ab  for  the  area  of  an  ellipse,  the  elemental  area  dA  of  the  skewed  ellipse 
on  axes  At^  ,  At^'  is  determined  by 


dA  -  •atN,4V-'i4V-d4V>'4V-i4Vi 

dA  -  7T  (A  t  •  dA  t  •  +  At  *  dA  t  '  _  dAt  1  dAt  ») 
r  m  M  n  n  m 


iG-14) 


To  reduce  the  niegration  from  a  double  to  a  single  summation,  we  evaluate,  in  effect. 


a  = 


/Sir  a  i 


//  exp- 


r  *» 

r 


(At  •)* 
m  ' 


3  a 


dAt  *  dAt  f 
n  m 


(G-15) 


for 


°<‘V<Ay 

0<AV<A<„' 


(G-16) 


The  exponent  of  e  is  the  constant  value  of  any  point  on  the  elliptical  ring  dA  as 
dAtn'  and  dAtm'  approach  zero.  But  the  total  area,  and  its  associated  joint  density 
function,  can  be  considered  a  concentric  series  of  such  rings;  as  the  maxima  At^' 
and  AtM'  move  from  zero  to  A  and  B.  For  example,  we  have  the  relationship 


2Tf"«  "  I  <G-17> 

M 

This  permits  us  to  generalize  Equation  (G-14) 

dA  - » 4  AvdA\  +!av  dav  -dAv  dAtm'>  ig-i8> 

But  from  Equation  (G-13) 


_1_ 

vsr 
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giving 


/ 

dA  =  7T  K'TAt..'  dAt  !  +  —  At..1  dAt  '  -  dAt  '  dAt 
V  M  m  /g  N  n  n  m 


) 


(G-19) 


The  integral  has  now  become 


a 


f  exp 


■fZTto  A 
x 


r  ^v>2 

3a2  a2 

X  X 


dA 


a  = 


V"3  V  O  2  A 

X 


/  exp  - 


r,A.n')2  (A«m') 

sa  2  a  2  ' 
Lx  x 


Let 


x  i/1  At.  '  dAt  '  +  ~  At*  dAt  ’  -  dAt  '  dAt  '  J 
y  M  m  N  m  m  n  J 


(At  f)2  (At  ')2 
'  n  ,  m  7 


=  S 


30 


(G-20) 


(G-21) 


be  a  general  ellipse  whose  value  is  constant  for  any  given  ellipse.  Since  it  will  not 
matter  whereon  any  given  ellipse  an  evaluation  is  performed,  i.e.,  At  1  =  At  ' 

At  '  -  0,  and  At  '  =  0,  Atm'  =  At^'  are  valid  maximum  points  on  the  ellipses  and 
we  may  write 


a  =  — 


VTo 


x 


r 

r  21 

(Ayr 

J  exP 

2 

30 

X 

At  1 

x  — —  dAtn'  +  J/3exp 
v~3 


r  ,AtM,»2Max"1 


AtM'  dAt 


m 


-f  exp[-S2]d^tu'  dAtm' 
The  last  term  may  also  be  written 

f  exp 


(G-22) 


\  <*M'>2' 


dAt  ’  dAt  » 
m  n 


It  is  obvious  that  the  expression  in  the  bracket  is  a  finite  integral.  Hence  in  the  limit 

,2  “ 


lim 

dAy  *■  0 


/exp  - 


dAt' 

m 


dAt  1  =  0 
n 


(G-23) 
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% 

& 

as- 

M 

? 


% 


I 

* 


Is 


the  evaluation  of  a  proceeds  an  follows 


1  r  3 

a  - -  —  exp 

75  o  2  75 


<a‘n’> 


1  y &  <AtMrM«  A  Atm' 

—  -5-  exp - s -  d  1 - 5" 

*  0  2  o*  j  V  “x 


The  limits  a'  and  b'  of  At^'  and  At^'  are  discussed  later.  Then 


1  (  Max 

“  =  -*eXP[—^ 

\  x  . 


-2  exp  I- 


Max 


(G-24) 


(G-25) 


c  -  +  ±|Lxp 


(G-26) 


It  is  interesting  to  note  that  for  special  cases  of  a  zero  ellipse,  and  an  all-inclusive 
Infinite  ellipse,  the  limits  and  probability  would  be 

a*  =  b'  =  0,  a  =  0 
a>  =  b’  =  of  =  1 

which  checks  our  knowledge  of  boundary  conditions. 

Since,  as  has  been  pointed  out  for  the  joint  probability  density  functions  of  this 
system, 

(a')2  =  3(b')2 


we  have  finally 


a  =  1  -  exp 


2 

=  1  -  exp  -  2 

(G-27) 

a  6 

t-  x  J 

3(7 

L  x  -» 
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or,  more  generally,  for  any  point  on  the  general  ellipse  S  where 


32  ,  ^q'*2  ,  <Atm')2  ^ 


At  At  At  At 
n  _ m  _ m _ n 

3/2  a  2  3/2  a  2  3/2  a  2 

X  X  X 


(G-28) 


a  =  l-e' 


and  the  points  a'  and  b'  are  the  special  "max"  values  of  At.,'  and  At-.'  which  give 
the  values  of  S. 

When  S  is  unit}'  we  have  the  natural  standard  deviation  ellipse.  The  probability 
that  the  position  fix  lies  within  or  on  this  ellipse  is  then 

o  =  l-e"1  =  63.22% 

For  concentric  ellipses  for  other  values  of  S,  the  elliptical  equation  can  be  considered 

9-8  n  r,  n  n 


wnere 


_  (Ay-)2  (Ay)2  At/  Atm2  Atm  Atp 

S2  3  a  2  s2  a  2  S2  3/2  a  2  S2  3/2  a  2  s2  3/2  a  2 

X  X  X  X  X 


At  "  =  SAt  ’ 
n  n 


(G-29) 


At  "=  SAt  • 
m  m 


and  the  unprimed  Atn,  Atm  are  also  changed  by  the  factor  S  .  Thus  we  may  consider 
any  size  S  ellipse  we  wish. 

To  relate  this  situation  to  the  space  domain  it  is  recalled  that 

n  =  At  r 
n  n 

m  =  At  r 
m  m 

a  -  T  o 

n  n  x 

o  -  r  a 

m  m  x 

Utilizing  these  quantities  in  Equation  (G-29)  gives  the  result 


-  S 


3/2  a  3/2  a  3/2  0  a 

n  m  m  n 


(G-30) 
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The  values  of  a  and  b  widch  were  used  in  Appendix  I  to  derive  the  transforms  are 
then  given  by 

a  -  s  rm  o 

n 

b  =  s  -im  o 

m 


These  are  pseudo-  or  orthogonal -elliptical  quantities  which  must  be  transformed  in 
accordance  with  Appendix  I  to  get  the  true  spatial  picture. 


APPENDIX  H 


DERIVATION  OF  THE  TRUE  ELLIPSE  TRANSFORMS  FOR  FOUR  OBSERVERS 

The  general  problem  is  to  derive  an  expression  for  the  true  ellipse  in  terms  of 
its  oblique  axes  such  that  a  complete  picture  of  the  spatial  probability  contours  is 
obtained. 

It  is  shown  in  Appendix  A  that  the  valr?  of  the  distance  R  from  a  fix  point  P  to  any 
displaced  point  ?'  is  given  by 

R44  =  — (m^  +  n^  -  2  rnn  cos  9  )  (H-l) 

sin  9  x 

x 


where 

6  -  crossing  angle  of  the  LOP's 

m  -  error  or  displacement  of  the  LOP  for  the  station  pair  AB 
measured  along  the  line  m  at  right  angle  to  rm  . 

n  =  error  displacement  of  the  LOP  for  the  statior  pair  BC 
measured  along  the  line  n  at  right  angles  to  rP  . 


It  is  also  shown  in  Appendix  F,  Equation  (F-26),  that  there  exists  an  elliptical 
mathematical  relation  of  a  constant  contour  of  probability  such  that 


The  question  is:  'Wheat  the  values  of  n  and  m  which  satisfy  this  equation  are 
substituted  into  equation  H-l,  can  we  obtain  the  "max"  and  "min"  of  the  resulting  R, 
and  what  are  these  values?" 


To  answer  the  first  part  of  the  question,  the  variations  of  R  are  studied  in  an 
example. 


First  m  is  replaced  in  Equation  (H-2)  by 


(H-3) 


(H-4) 
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R  =  1ET~  ^ q2(1  "  b2//&2)  +  fa2  -  2cn  -  n2/a^  cob  0  (H-5) 

For  the  example  chosen 

30° 

8 
2 


and 


R  =  J  3n2  -  9.  8n  /T^iT/S  +  8 

The  results  of  this  calculation  are  shown  in  Figure  H-l.  This  demonstrates  that 
there  is  a  value  **f  R  =  f(n)  which  in  the  usual  orthogonal  sense  will  contain  a  maximum 
and  Tnininr’m. 


f? 


Figure  H-l.  A  Plot  of  an  Example  of  Calculating  R 
to  a  Point  on  the  (Oblique)  Ellipse 
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« 


To  derive  the  general  expression  for  all  maximum  and  minimum,  we  return  to 
liquation  (H-4).  By  inapt  ction  of  the  ellipge  of  Equation  (H-2)  projected  onto  the 
oblique  system  n,  m  (as  n  Figure  6)  it  is  seen  that  Rmax  ^  Rmin  obtained  when 
dR  =  0.  But  since  dn  is  independent  of  dR,  if  dR  is  zero 


dR 

*ST" 


0 


From  equation  (H-5), 


-i/2 


2n  -  2— ^  n  -  2b 

SL 


2 

Equating  to  zero  and  solving  for  n 


proceeds  as  follows 


A 

7/ 

.  b  2 
+  n  cos 

sT 

°x  '  b(*  "  008 

1 

|c r 

+  2  -y  cos 

6  n2  -  b  cos  9=0 

a  / 

a 

=  b  cos  0 

0 

-  2  b,-  a  cos  & 

u 

a  / 

=  0 
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Lquaring  both  sides 


2  n2W.  b2\  .2  2  .  b2  2  a  2-  ,  b2  2  4 

b  ^1  — ^ j  ^1  — g y  b  cos  0  -  4  -£  cos  en  +  4  cos  0  n 


2  4  .  b  2  2  ±  ,  2  2  a 

cos  8n  -4-^  cos  0  n  +  b  cos  0 

a2 


2 

/,  fc2\  n4  /.  b2^  ,  b2  2  4  .  b2  2  2  .  ,  2  2  a 

U  - -j  J  — ^  U  -  “Jj  =  ^  -j  cos  9n  -  4  -j  cos  8  n  +  b  cos  8 

\  a  /  a  \  a /  a  & 


B)-‘- ('-9 


2  i  b  2  a  2  ,  1 2  2  a  A 

n  -  4  cos  8  q  +  b  cos  0=0 


°“2<,+  (1'?)  ]?'  . 

?[4ScM2<,+(i'S)]a4* 


,  b2  2  n 

4  cos  8  + 


1 4  ^  cc©2  8  + 
a2 


(■  ■  S)  ]■ 

(-Sft 


2  +  b2  cos2  0  =  0 


2  +  b2  cos2  0=0 


Solving  for  n  using  the  binomial  theorem 


,2  2  u 
b  cos  0 


14  2  A 

_  n  .  n  +  — - 

4  cos2  0  + 
a*5 


r 

i  ±i/ 1  - 


. .2,  2  2 

4  b  /a  coi  0 


n2  = 


4b2/a2  cos2  0  +  (1  -  b2/a2)2 
2/a2 


n2 

~2 

a 


▼  4  b2/a2  cos2  0 


(H-6) 
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For  computation  purposes,  let 


1  -  b2/a2  _  a/b  -  h/p 

2  b/a  cos  6  2  cos  & 


(H-7) 


Then 


but  from  Equation  (H-2) 


and 


(H-8) 


(H-9) 


From  Equations  (H-8  and  H-9):  the  values  of  m  and  n  which  give  Rmax  and  Rmin 
are 


n  =  ±~Jl±Jl  -  -HI  (H-10) 

V2  "  *  1  +x“ 

m  =  ±—  Jl+Jl-  — 1— ^  (H-ll) 

/2  "  f  1  +x2 

Insertion  of  these  quantities  into  Equation  (H-l)  presents  a  problem  in  the 
selection  of  polarity  for  the  term  mn  cos  6X,  This  has  been  solved  by  graphical 
analysis  which  shows  that  in  the  determination  of  Rmax.  m  and  n  are  always  of  opposite 
polarity  regardless  of  whether  the  pseudo-ellipse  is  plotted  with  its  jcifljor  «xis  «icng ' 
m  or  n.  Similarly,  graphical  analysis  shows  that  for  the  point  on  the  true  ellipse 
which  represents  Rmin  m  and  n  are  of  the  same  polarity. 
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tm  -.m.  *> 

re/<f  •  ,••?■*■ 


Further  analysis  with  regard  to  the  selection  of  polarity  shows  that  when  a  >  b, 
Rmax  is  obtained  from  the  combination  of 


°2  =  K1  V77?.  x2) 


Conversely  when  b  >  a,  Rmax  is  obtained  from  the  combination  of 


m2  _  1  ( 
2  ' 


+*  1  - 


\  T  1- x/ 


Thu«  the  formulas  for  the  true  major  and  minor  axes  become: 


For  a  >  b 


2  sin  d 


^fr5)+b2tf-7b) 


+  2  ab+l  (l  ±*/l 


Simplification  of  the  last  term  reduces  to 


■rbH^rSH 


a2  -  1  ,2 

8  max  ”  2  »in2  0 


^  +b*  (‘  -V*  -it?) 


2  ab  cos  6 


(fl-12) 


Similarly,  for  a  >  b 


min  r>  ,  2  ^ 
2  sin  6 


J-V77*)  *“■(!•  V77^ 


2ab  cos  0 


ffi-13^ 


and  for  b  >  a 

R2  =  - V- 

max  2  siiT0 


4 b2  (  'iFtt 


2/ 

X/ 


2ab  cos  ft 

/777 


(H-14) 


R 


min  „  .  2  a 
2  sin  8 


a2  (l  +Jl - O  -'•b2  [l  /l  -  — O', 

\  t  1  +  x  /  V  ' 


1+J) 


2ab  cos  9 


/ 1+7 


(H-15) 


If  the  angle  of  shift  of  the  pseudo -major  axis  is  4'  which  is  measured  from  that 
axis,  m  or  n,  to  the  true  major  x~axis,  this  angle  can  be  calculated  from  the  above 
information.  Basically  we  are  dealing  with:  for  a  >  b. 


The  n  value  of  R 


COS  4/  = 


max 


(H-16) 


max 


2 

cos  y 


aV2  ( 77) 


2  sin  8 

which  reduces  to 


cos  9 

2 


sin2  6 


cos  i|/  = 


v1  V^TT?) 


1  +JTTJ—  +  b2/s  (l  -ind  +  .2.i>/a.cos  * 
1  1+*  \  *  /777 


(H— 17) 
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mmmmmirj*-  ■  mam**™-*** 


Similarly  for  b  >  a 


sift"’  6 


2  . 
COS  4/ 


1  +  3?  > 


.V  iJi 


-T7?-)+(1+lFr7?)  +  ^ 


2a/b  cos  0 


(H-'it) 


Briefly,  regarding  symmetry,  if  we  had  started  the  original  solution  in  terms  of 


dB  _ .  n 
0 


the  resulting  solution  would  have  given 


1  ,  /  .  4  a2/b2  cos 2  0  \ 

r*W 


Then  let 


1  -  a2/b2  _  b/a  -  a/b 

2  a/b  cos  0  2  cos  0 


which  would  become 


Ui  —  f -  . 


From  this 


but  we  also  have  the  dR/dn  =  0  derivation  that 


Therefore 


Continuing  this  derivation  gives 

x  =  ±y 

which  for  x  =  -y  substantiates  the  condition  obtained  in  the  derivations 

x  - - l -  (g/b  -  b/a) 

2  ccs  &x 

y  =  rdnr  <f =  -Toknr  <|-b/a> 

X  X 

It  is  desirable  to  provide  transforms  for  the  major  and  minor  axes  such  that, 
given  values  a,  fc,  and  8%,  the  true  values  of  the  major  aad  minor  axes  and  the 
orientation  of  these  axes  can  be  determined  directly  and  simply.  Such  a  transform 
would  be  the  ratio  of  values  of  the  true  major  axis  to  the  pseudomajor  axis  a  or  b. 
Since  the  original  maximum  axis  value  could  bo  a  or  b,  wo  will  call  this  and 

the  true  value,  R'msx*  The  transform  (squared)  then  becomes 

for  a  >  b 


2b/ a  cos  6 

X 


(H* 


m 


It  can  be  seen  by  inspection  of  these  formulas  that  if  the  b/a  of  Equation  (H-19) 
equals  the  a/b  of  Equation  (H-20),  then  the  x  of  Equation  (H-20)  will  be  the  negative 
of  the  x  of  Equation  (H-19),  but  the  squares  of  these  two  will  be  equal 

*"XH-20*  =  (+  XH-19) 

Hence  the  transform  functions  are  equal  for  given  values  of  a/b,  b/a,  and  It  is 
most  convenient  in  visualizing  these  transformations  to  picture  the  pseudoellipse  as 
having  its  major  axis  coincident  with  the  corresponding  axis  in  the  oblique  system  m, 
a.  Thus  if  a  >  b,  we  would  draw  the  ellipse  to  be  transformed  with  n  as  the  major 
axis.  Likewise  tf  b  >  a,  m  would  be  considered  the  major  axis  of  the  pseudoellipse. 
3y  similar  development  the  square  of  the  minor  axis  transforms  becomes 


“U  1 


— 4-  i  4  (i  -J'  ■~14) +  f 

2  sin'3  B  j  b2  \  >  I  +  x2/  \ 


+J 1  - 


2 a/b  cos  0 1 

41  +  x^ 


(H -21) 


and  for  b  >  a 


CR»  .  ) 
min' 


2  sin2  6 


1+V7T77+S(-Vr:TT 


(H-22) 


2b/a  cos  d 


1  +x 


ai  d  since  the  same  conditions  prevail  with  respect  to  the  x  values,  the  symmetrical 
^quality  of  these  two  equations  is  also  established. 

Returning  to  the  determination  of  ,  it  can  be  seen  by  comparing  Equations 
(H-17)  and  (H-19)  that 


And  finally, 


cos  \j/  = 


1  +*/ 1  - 


cos  \f/  = 


2  (Major  Axis  XFrm) 


(H-23) 


(H-24) 
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APPENDIX  I 

DERIVATION  OF  THE  TRUE  ELLIPSE  TRANSFORMS  FOR  THREE  OBSERVERS 


As  in  the  case  of  four  observers,  this  derivation  starts  with  the  value  of  R  as 
given  by  equation  (F-26) 

R2  =  — ^ - (m2  +  n2  -  2mn  cos  6^  (I-l) 

sin  8 

x 

It  is  shown  in  Appendix  J  that  there  exists  an  elliptical,  mathematical  relationship  of 
a  constant  contour  of  probability 

2  2 

5_  +  SL  _  ”95  =  i  n-2) 

2.2  ab  1  (1 

a  o 


that  is 


n2  =  -~r  (m  i  i/4b2  -  3m2 


I 


giving,  by  substitution, 


sin  8  |_ 


m2  +-^U  i^b2  -  3m2^  -  m|  ± ^4b2  -  3m2 


(1-4) 


where 


C  =  cos  9 


H  '  *  ^2  (m  *  -  3t"2) 


V4b2  -  3m^C  (M 


Again  proceeding  as  in  the  case  of  independent  observers,  Rmax  and  ^min  sre 
obtained  when  dR  =  0.  Since  dn  or  dm  is  independent  of  dR,  we  are  also  solving  for 


dR  1  „  .  a  -1/2 
to  =  2  ®  8,n 


[2ra+-4 

L  4bZ 


x  2(ir.  +  y4b  -  3m  )  (l  - 


.  ■  ■ ...  ■  IIWU.„,  ,  l 

274b2  -  3m 2) 


i-2 


Sohrlnf  for  m  proceeds  aa  follows. 

To  reduce  repetltlri,  let  x  =  v4b2-  Sm2 
Then,  from  equation  (1-6), 


2m  +  — 
2b 


2ns  + 


(*  +  -|(m  +  x)C.C 


2mx  +  (x  -  3m)  1  (m  + 
L2b 


2mx  +  ^(mi  ♦  ,2,  -  i  c  rax  .  Ill  (m2  +  m)  +  Sac  m2  .  a  + 


x)-mfc]-| 


C(mx  +  x^)  =  0 


C(mx 


2mx+--mx  +  i  x2-£cmx-|-^r  m2  Jl 

ou*  b  2 . 2  9 


u  «  A  o  a  a  9 

2^2  ”  -2^mx  +  3bCm 


.iPmv  a  ^  2 
-  •g  C  mx  -rCx  =0 


x(2m^m-fcr,-!^m-fcm)+x2^-|c^  + 

ra2^C'^)'° 

™(2'^c'A+x2^'bc)+m2(3fc-!^)-» 

m(2-2|c-^+^.ic).3m2(Zc.ac^o 


xm  2  -  2  ~  C 


4b4  -  12b2v  +  9v2  =  4k2b2v  -  3k2>  2 


3(3  +  k2)v2  -  4b2 (3  +  k2)v  +  4b4  =  0 


3v2  -  4b2v  +  =  0 

3  +  k2 


4b2  ±J  16b4  -  48  — — — - 
T  3  +  k2 


v  =  m  = 


-•■n-ra 


(I-H) 


Finally  to  get  this  into  a  form  comparable  to  the  four-observer  case,  let 


A  _  JL 

5  3 


giving 


2  2.2 


/  2 
19  o  p 

V  15  V. 

(‘•Ftt?. 


(1-15) 


d-16) 


as  the  coordinates  of  points  on  the  true  ellipse  which  are  the  end  points  of  the  major 
and  minor  axes . 

'‘'he  n  coordinates  are  obtained  directly  from  considerations  of  symmetry.  If 
initial  substitution  in  equation  (1-6)  had  been  made  lor  m  instead  of  n,  and  the  equation 

dn 

solved  for  n,  the  resulting  equations  are  the  same  as  those  above  if,  also,  the  following 
replacements  are  made 

m  — »  n 
r. --*>m 
a  — >  b 
b  — ►a 


Then 


2.2  2 A  j.  1  ' 

3  3  \  ’  1  +  J2> 


a-n) 


where 


t2  V  »  a2> 

•  /  O  \  O 


(1-18) 


For  a  quick  demonstration  of  the  correctness  of  these  equations,  consider  the  special 
case  of 

6  =  90°  (i .  e . ,  orthogonal  axes) 
a  =  b 


j.2  =  (2  -  0  -  1)  =  1 

3(1  -  0)2  3 

2  2 ,2n  .  1, 
m  ~  3  b 


m  =  ± —  and  ±b 
V3 


n  =  ±  —  and  ±a 
V3 


Remembering  that  when  a  =  b, 


2  2 

m  n  mn  _  . 

72  +^"Tb  ' 
b  a 


is  a  skewed  ellipse  of  45°  given  by 


£sai+aaf  =  i 

c  2 

lb')  (a'j 


(1-19) 


wt^ere 


<b’)2  = 


cos2  45°  ,  cos2  45°  sin  45°  cos  45° 


(1-20) 


X7X7X  ' 21,2 

2b2  2b2  2b2 

b'  =  -h  b 


(1-21) 


«-.»2  1  2.2 

<a  )  -X7X7X  sb 

2a2  2a2  2a2 


a-22) 


a'  =^a 


a-23) 


These  values  are  shown  on  the  following  sketch  which  illustrates  and  verifies  the 
special  case. 


Another,  more  typical,  example  was  worked  out  analytically  and  then  checked  with 
graphical  solution.'  The  parameters  selected  were 


a  =  0.25 
n 

a  =0.6 
m 

s  =  1 
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giving 


*  »#' ».  -  °-306 


»=9l|f0m  = 

We  then  wish  to  find  the  locus  of  the  points 


0.735 


m2  ,  _»L _ -55S -  =  1 

(0.7 35) 2  (0 . 306) 2  (0.735  *  0.306) 


onto  the  oblique  anon  to,  n,  »  The  «Ud ^solution  IB  shown  In  Figure  1-1 . 
plotted  values  were  obtained  from  the  relations  p 

n  =  ^  (m  ±1/4?  -  3n?) 

The  coordinates  from  this  graph  for  the  terminus  of  Rmln  and  Rmax  ar® 

forRmfl:  m  -  0.780,  n  =  0.043 
max 

forRmin:  m  =  0.325,  n  =  0.348 
To  solve  for  the  coordinates  analytically 

—  =  0.416,  cos  20°  =  0.94 
b 

_  2  -  2X  0,416  X  0.94  -  0.173.  =  _  Aill 
^  ~  Vs  (0.173  -  2  X  0.416  X  0.94)  V3 

|2  =  0.975 

m2  =  §  b2  (j  *  1^2<1  *  °-703) 

m  *  1.063b,  0.445b 
m  *  0.782,  0.327 


Likewise,  for  n 


2  -  2X  2.4  x  0.94  -_jL76_  _  _3  82 
/3  (5.76  -  2  X  2.4  X  0.94) 


/  ? 


t2  =  14.6 


n  = 


;V 


2„  7 


6 )  3 

n  =  0.146a,  1.145a 
n  =  0.0447,  0.350 


x  0.306  (1  +  0.968) 


And  the  coordinates  become 


R  :  m  =  0.782,  n  =  0.0447 
max 

R  ,  :  m  ~  0.327,  n  =  0.350 

in  * 


Therefore,  we  are  dealing  with  something  less  than  2%  error,  including  graphical 
and  analytical  calculations,  and  the  formulas  are  verified  for  this  example. 


Providing  transforms  for  the  true  values  of  the  major  and  minor  axes  is  more 
complicated  than  for  the  previous  case.  First,  the  quantities  a,  b  no  longer  repxesent 
the  major  and  minor  axes  of  the  pseudo-orthogonal  ellipse.  They  still,  however, 
represent  convenient  quantities  for  determining  the  desired  transforms. 

Returning  to  the  general  form  of  the  radial  vector  to  the  true  ellipse,  and  substi¬ 
tuting  the  derived  values  for  the  true  coordinates  m,  n,  givps 


R2  = 


sin2  & 


2  ,  2 
3b 


itflW tS)-' ■ 


fl-24) 


Defining  the  transforms 


XFM«=f  P"25> 

XfMH  =F  (a>b)  “7(b>a)  (1-26) 

The  neede  i  relationships  R/a  and  R/b  are  determined  from 


Since  the  key  factors  to  those  solutions  are  the  quantities 


and  since  one  abscissa  can  provide  either  a/b  of  |  or  b/a  of  i ,  Figure  1-2  is  a 
representation  of  either  radical .  The  selection  of  the  parameter  6X  gives  samples 
through  the  crossing  angle  spectrum. 

The  maxima  (of  unity  value)  ars  obtained  for 

b/a  or  a/b  =  2C  (2  cos  0  1 

The  minima  (of  zero  value)  are  obtained  for 

b/a  or  a/b  =  -C  +V?  +  2 

The  case  for  6  =  60°  is  very  interesting  and  perhaps  eveu  anomalous,  The  full  sig¬ 
nificance  of  this  curve  has  not  been  explored  for  this  document  exceot  that  it  passed 
ad  validity  or  application  tests  developed  during  this  writing. 

The  next  task  is  that  of  making  the  various  polarity  decisions  associated  with  the 
radicals.  There  are  two  distinct  decisions  to  be  made:  (1)  the  polarity  of  the  radical 
in  determining  the  value  of  the  coordinates,  and  (2)  the  relative  valus  of  the  coordinates 
themselves  with  respect  to  each  other.  This  is  a  much  more  difficult  task  than  for  the 
independent  case  of  only  one  radical  value  per  a.  b,  c. 

The  tools  ard  criteria  for  making  these  decisions  were 

a.  Figure  1-3 


Figure  1-3 .  Geometrical  Polarity  Decision  Criteria 


■fern 


b. 


2  2 

m  n_  mn  _  . 

2  ,  2  '  ab  1 

a  b 


c .  Some  trial  and  error 


The  results  are  tabulated  in  Table  1-1  which  seemed  to  be  about  the  best,  though 
by  no  means  the  only,  way  of  depicting  these  polarity  decisions .  Since  the  quantities 


2C  and  -C  + 


fc2  +  2 


can  both  be  greater  or  less  than  one,  ail  possible  cases  are  charted.  The  idea  used 
in  deriving  the  chart  was  to  find  a  starting  point  for  which  all  the  factors  wore  known 
and  then  reason  what  must  happen  from  then  on.  (The  author  is  sure  there  is  a  ~urer 
mathematical  or  classical  way  of  doing  this,  but  the  idea  was  successful.)  The  starting 
point  used  is  b  =  a.  Examination  of  the  elliptical  equation  shows  that  the  choice  of 
radical  polarities  is  limited  to  the  extent  that  they  must  be  the  same  for  m  and  n  at 
b  =  a.  The  next  significant  factor,  obtained  by  examining  Table  1-1,  is  that  Bx  =  60° 
is  a  dividing  point  for  these  polarities  and  the  cases  of  9X  >  60°  must  be  examined 
separately  from  0X  <  60° .  Finally,  the  separation  of  m  and  n  for  XFMax  and  XFMin 
is  obtained  by  almost  any  test  points  within  the  regions  specified. 

Having  thus  obtained  a  complete  initial  set  of  polarities  for  the  radicals  at  b  =  a, 
the  rest  is  the  intelligent  guess  that  as 


1  +  or  t2) 


goes  through  zero  at 


the  polarity  associated  with  the  radical  must  reverse  as  a/b  and  b/a  are  varied 
through  their  interlocked  values . 

The  significant  factor  in  determining  the  polarities  of  m  and  n  relative  to  each 
other  is  the  value  2C,  for  at  a/b  or  b/a  =  2C,  m  or  n  are,  respectively,  zero.  This 
would  be  intuitively  the  point  at  which  a  coordinate's  polarity  would  change.  Verifica¬ 
tion  of  this  change,  and  determination  of  the  manner  in  which  the  polarities  changed, 
were  accomplished  by  graphical  analysis,  with  the  results  shown  in  Table  1-2.  It 
should  be  remembered  that  as  far  as  the  transform  formulas  are  concerned,  we  need 
know  only  whether  the  relative  polarities  are  the  same  or  opposite. 

Determination  of  the  angle  of  major  axis  shift  ^  is  performed  in  a  manner  similar 
to  the  analysis  for  four  observers.  It  proceeds  as  follows 
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For  b  >  a 


or 


2  [value  of  m(for 
co«  j - 

RMax 


Max 


where  the  polarity  of  the  radical  is  determined  in  the  same  manner  (using  Table  1-1) 
as  for  the  XFMax*  The  angle  V  would  then  be  measured  from  the  m  axis . 


Similarly,  for  a  >  b 


CO  8  \p 


and  4/  would  be  measured  from  the  n  axis . 
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TABLE  I-i 


POLARITY  DECISION  CRITERIA  FOR  THE  RADICAL 
DETERMINATION  OF  RADICAL  POLARITY 

(*  =  -C  +  V? "+2  ) 
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APPENDIX  J 


DERIVATION  OF  THE  JOINT  PROBABILITY  DENSITY  FUNCTION 
FOR  THREE  OBSERVERS 


The  basic  problem  is  to  find  the  true  degree  of  dependency  or  correlation  between 
the  time  difference  pairs  which  have  one  observer  in  common.  Suppose  observer  B 
were  the  center  or  common  observer.  Further,  for  notation,  let  the  errors  or  devia¬ 
tions  from  true  time  readings  be 


TA-‘A  =  ilA 

(J-D 

TB  ”  *B  =  **B 

(J-2) 

Tc  -  -q  -  At., 

(J-3) 

where  T^,  Tjj  and  Tq  are  the  true  instantaneous  values  of  time  of  arrival  and  t^, 
tjj,  t£  are  the  observed  instantaneous  values. 

Figure  J-l  presents  one  example  of  a  possible  time  configuration  with  probability 
density  functions  sketched  for  each  observer.  The  probability  that  A  has  made  an 
error  of  magnitude  between  At^  and  AT^-t-dAty  is 


Similarly, 


p(At  )  » - —  exp 

aAV2» 


AT, 


20. 

l  A 


dAt, 


P(Aty)  = - — •  exp 

(Tg  V2T7 


r  2' 
AT_ 

C 

2cb2 


dAtB 


P(Atr)  = - __  exp 

0Q  V27T 


AT, 


2  cr 


C  J 


dAt, 


Defining  the  quantities  At^  and  At^  as 

AtN  ^A  ”  TB)  “  (tA  "  V 
AtM  =  "  TC>  "  ^ 


(J-4) 


/T_c' 


(J-6) 


(J-7) 

(J-8) 


the  goal  of  this  derivation  is  to  determine  the  joint  probability  function  f(At^,  At^). 
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Figure  J-l.  An  Example  of  Probability  Distributions  for  Three  Observers 
Subetilut  ng  from  equations  (J-l)  and  (J-2) 


and,  finally. 


sf 

£ 

11 

«T 

(J-9) 

*m* 

(J-10) 

*A-*N  +  AtB 

(J-ll) 

Atc  "  alb  '  a*m 

(J-l 2) 

Considerable  confusion  and  discussion  can  arise  on  the  relative  polarities  of  these 
quantities  and,  as  indicated  in  Figure  J-l,  the  possibilities  include  all  combinations 
of  polarities  o  1  the  three  quantities.  Since  we  are  concerned  only  with  AtA2  and  At^2, 
thr  se  combinations  can  all  be  reduced  to  two:  Atg  is  either  aiding  or  opposing  the 
<;  entities  AtN,  AtM. 

The  next  important  probabilistic  concept  is,  that  to  obtain  a  specific  Atjj  and  At^. 
we  must  also  obtain  specific  values  for  AtA,  Atg,  and  At^.  The  probability  of  doing 
this  is  the  probability  of  obtaining  these  three  errors  simultaneously,  and 


p(AtN’  AtM>  =  p(AtA)  p(AtB)  p(AtC) 


(J-l  3) 


where  Au ,  Atg  and  At^  represent  those  quantities  necessary  to  obtain  a  desired 


AtN  and  At 


-M- 


Substituting  from  equations  (J-4),  (J-5),  and  (J-6) 


D(AtN’  AtM> 


3/2  6XP 

Wc<2n> 


/At, 


V 


At, 


2a. 


> 


2aB2  2SC2/J 


dAtAdAt£dAtc  (J-14) 
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It  is  assumed  at  this  point  that  characteristics  cf  the  three  observers  are  so  nearlj' 
alike  that 


*A  =  aB  =  ffC  *  * 


P(AtN’  AtM*  a3(2ir)3/2  6Xp 


(At  2  +  At^  +  At  2>] 

” ~  "17 -  J  dAtAdAtEdAtC 


(J-15) 


Now,  substituting  for  At^  and  Atc, 


p<AtN' Atw) "  7< W71  eTP 


K2  +  <AtN4V  +  <AtMiAtB>:l 


jdAtAdAtBdA,C 

(J-18) 


To  obtain  the  total  probability  for  the  occurrence  of  a  given  AtN  and  Atj^,  we  must 
add  the  probabilities  for  all  the  mutually  exclusive  collectively  exhaustive  ways  in 
which  these  AtN,  Atj^  can  be  obtained.  This  is  done  by  allowing  Atg  to  successively 
take  on  all  possible  values  from  -« to  +oo. 


dAtAdAtC  r°° 

p(At  ,  At  }  =  — - 3 7=-  J  exp  -  - 

aJ( 2ty/l  L 


O  (2tf) 


V  4  <AtN  4  AV2  +  <AtM  4  AtB>2l 


jdAtB 


(J- 17) 


Since  At^,  Aig  and  At^  are  independent  and  since,  for  example,  At^  and  Atg  are 
kept  separated  by  AtN, 


—^AtN  \+- 
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and  we  can  conclude 


dAt, 

dAt^  = 


dAiA  =  dAtg  =  dAtc  =  dAt 


Equation  (J-l?)  can  be  rewritten  as 

p'Atjj,  Atjj)  _  j  90 

JZZJ72 


*B~  +  <AtN  *  V 


20* 


(dAt)  cT  (2  ST) 

If  we  then  define  the  joint  probability  density  function  as 

pf'V  -V 


dAtg  (t>-18) 


f<AV  A^)  - 


(dAt)' 


( J- 1 9) 


the  result  is 


00 


f(<V  AtM’  =  (,3^3/2  /„  exp  - 


atB  +  (4tN  *  AtB*  *  (AtM  *  AtB*  1  .. 

- —2 - J  ditj,  (J-20) 


For  those  interested  in  a  more  classical  verification  of  this  relationship,  the  problem 
is  one  of 


given  P(x,  y,  z) 

s 

"\  find  P'(u,  v) 

where  x,  y,  and  z  are  independent  variables  and  u  and  v  are  some  functions  as 

u  =  u(x,  y,  z)  and  v  =  v(x,  y,  z)  (J-21) 

For  any  particular  u,  v  there  is  a  curve  in  x,  y,  z  space.  The  probability  of  being 
in  the  incremental  area  (u  +  du  -  u)  by  (v  +  dv  -  v)  is  also  the  probability  of  being  in 
the  volume  x,  y,  z  which  is  formed  by  the  locus  cf  u(x,  y,  z)  and  v(x,  y,  z)  for  each 
given  u,  v.  If  the  curve  is  single-valued  with  respect  to  the  z  axis, for  example,  we 
can  integrate  over  z  to  find  the  total  probability  of  being  in  the  volume .  To  find  these 
curves  as  a  function  of  z,  the  equation  (J-21)  must  be  solved  to  give 

L 

x  -  x(u,  v*  s)  (J-22) 

y  *  y(u,  v.  z)  (J-23) 
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Then 


P'dudv  =  Y  ?[x{u,  v,  a),  y(u,  v,  z),  z]dA 
All  z 


("24) 


where,  using  vectc"  notation, 
dA  = 


^ + 1?  x  (l? \ + ^  | dudv 


fJ-25) 


and  mx»  #*y»  vx»  v„  are  unit  vectors.  When  the  cross  product  is  performed,  this 
reduces  to 


Therefore 


P1  -  /P[x(u,  v,  z),  y(u,  v,  z),  *]||2|j  -|f  Jjjda 

where  P'  is  the  desired  probability  density  function. 

For  the  example  at  hand 


where 


Evaluating 


a  o/o 
C'(2TT)''/t’ 


2 


X 


+  y2  +  z 
20* 


x  =  AtN  vAtg 

y =  AtM  +  AtB 

z  =  Atg 

u  =  (x  -  z)  -  AtN  +  Atg  -  Atg  -  Atjj 
v  =  (y  -  z)  =  At^  +  Atg  -  Atg  --  At^ 


lx.  8<*K  +  AtB> 
9U 


(J-26) 
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gives 


and 


5x  8y  8x  8y  I 

r 


QC 


P'  =  "o— 1— o7n  f  sxp 
trnvy'2 


which  is  precisely  the  3ame  as  equation  (J-20) . 


(Aijj  a  Atg)2  +  (A^  *  Atg)2  +  Atg2 


20 


dAtg  (J-27) 


Continuing  with  the  solution  of  aquation  (J-20),  when  the  exponent  is  expanded  and 
collected,  it  can  be  written  in  the  form 


3 

202 


*1*“ r*  ~it~ 


from  the  integral  tables  which  state 
00  r 

/  exp  -  (Ax“  +  Bx  +  C)  dx  -  exp 


B  -  AC 


] 


(J-28) 


Substituting 


1  7T  ( 

r<AtN’AV2  3  AtN2  +  At;/1 

_ 1  _ 

J  o  exp 

1  3/202 

1 

i 

sj 

to 

So 

3/2CT2 

(J-29) 


when  the  coefficient  and  exponent  are  simplified 


t/M  a  -  1  \  2/3(AtN  *  AtM  ~  AtMAtN*  I 

(  N’  AtM)  "  ^  6XP  ['  202  J 


(J-30) 


and  using  the  standard  deviation  per  pair  of  stations,  as  noted  in  Appendix  E,  i„g., 


a  2  =  20^ 


f(AtN,  At^  =  -j-j;  exp 


2/3(AtN2  <  At^  -  Ayy 
2 

u  J 


(J-31) 


(J— 32) 


which  is  the  formula  discussed  in  the  main  text. 
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It  is  interesting  to  note  the  direct  comparison  between  this  result  end  the  result 
of  Gaussian  noise  adder  circuits.  If  we  have  three  Gaussian  noise  inputs,  x,  y,  z 
with  standard  deviation  <T  for  sach,  and  use  an  adder  (or  subtracter)  circuit  element 
for  x,  y  and  z,  the  situation  symbolically  looks  like 


The  density  function  for  this  condition,  given  in  numerous  texts,  is 


f(u,v)  = 


W2  -h-o' 


exp 


2  Vo  yi-p 

where  p  is  the  correlation  coefficient  obtained  from 

uv 


1 

2  -  2ct,  +  v2 

“  2 

« 

.  twW  . 

I 


P  ~ 


a  a 
u  v 


However 


and 


hence 


uv=(y-x)(y-z)=y  -  zy-yz  +  xz 
uv  —  y^ 

?  =  S+?  =  f, 

l!V  =  o2  \ 


It  can  ai  so  be  shown 


(J-33) 
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Fes ee 


P 


g2 

V2oV2g 


1 

2 


Aad,  for  thi*  vidua  o i  p. 


f{u,  v) 


2  2 
(u  v 


g2 


which,  again,  l*  identical  to  the  derived  formula. 


I_£5l 
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APPENDIX  K 


MODIFICATION  OF  ELLIPTICAL  TRANSFORMS  FOR  A  DIVERGING, 
OBLIQUE  COORDINATE  SYSTEM 


As  stated  la  the  main  text  the  true  coordinate  system  about  the  point  P  is  not  a 
gridwcrk  of  lines  parallel  to  the  LOP's  at  the  assumed  true  point  P.  Rather,  each 
set  of  points  m,  n  which  satisfy  the  error  conditions  At^  Atjj,  lie  on  the  hyperbolas 
associated  with  Atj^,  Atj.  so  that  the  angles  between  the  m  and  n  axes  and  the  hyper¬ 
bola®  are  right  angles . 

If  we  (see  Figure  K-l)  also  construct  the  intersections  of  ±m,  *n  in  the  usual 
manner  on  lines  parallel  to  the  LOP’r,  then  these  points  and  the  points  described 
above,  for  example  m1  and  n',  and  the  true  point  P  all  lie  in  a  straight  line.  It  Is 
then  desired  to  determine  the  change  or  adjustment  of  the  values,  called  R  (or  C)  in 
the  text,  into  the  values  d  and  a.  This  is  accomplished  with  engineering  approxima¬ 
tions  in  the  following  manner. 

In  addition  to  Figure  K-l,  consider  the  following  enlargements  of  the  regions 
8  and  8  . 

x  y 


/ 
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Then  by  trigonometry 


$x  =  A  +  B,  Qy  -  E  +  D 


Sin  A  =  -£  Sla  B  =  — 
-  c'  ■  c 


sin  A  _  g 
sin  B  h 


sin  £  ~  sin  x’  sin  A  sin  80 

h  _  a  h  _  c 
sin  D  -  sin  p  *  sin  B  sin  90 

sin  E  _  g  sin  x 
sin  D  "  h  sin  m 

,TTe  then  exercise  the  engineering  approximation  that 

■ifiJE  o.  i 

s  in  M 


(K-l) 


(K-2) 


(K-3) 


(K-4) 


sin  E  _  g 
sin  D  h 


(K-5) 


Further 


sin  A  “  A,  sin  B  “  B,  etc. 


Therefore 


A  M  g  E  «  £ 

B  h’  D  h 


(K-6) 


A  ■,= 

B  D 


(K-7) 


A  +B  M  E  +D 
B  D 


(K"8) 


B  D 

T  =  F 

x  y 


(K-9) 


however,  as  before, 


c  _ J5 — 

(K-15) 

sin  90  sin  D' 

d  5 

(K-16) 

sin  y  sin  B* 

Bin  D!  siny  „ 
sin  B! 

D' 

B' 

(K-l?) 

and 


Finally 


D’  +  E'  =  Q% 
A'  +  B'  «  9Z 


To  convert  this  to  other  system  parameters  considor 

+  x  +  p  +  <f>'  =  360[f  *  180  -  0j 


(K-18) 


(K-19) 


(K-20) 


(K-21) 


(K-22) 


or 


and 


&  ~  180  +  0  -  (x  +  p) 

y  x 

5  =  x  +  p  +  0  “  180 
z  * 


Further 


=  90  -  x 

2 


n 


=  90  -  p 


(K-2S) 


(K-24) 


(K-25) 


GMj* 


A0  +  A0 

- g - -  »  180  -  (x  +  n) 

/a*  +  AS  \ 

(  mg-  »  (R  +  M)-180 

Equations  (K.-14),  (K20)  and  (X-21)  can  be  rewritten 

«  »  +  ££e3> 

C  y  +  2  +  2 

- 

X  X 


(X-26) 


(K-27) 


(K-28) 


w 


x  «  _ m  ~  n 

x  2  2~ 


(X~29) 


O  x  ^  x  -■-  b 

d  9x  “  2 


a  _ m  _ n 

“  2  2 


Finally,  to  a  very  close  approximation, 


-J5  , 

2  r  y 

m  m 


2  r 


(K-30) 


and  (with  6  given  in  radians) 


Q  +M+M, 

d  x  rn  rm 

a  "  g  .  m:m 


x  r 


n  m 


(X-31) 


and  so  forth. 
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